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Abstract 

During  the  August-September  2003  Autonomous  Ocean  Sampling  Network-II  ex¬ 
periment,  the  Harvard  Ocean  Prediction  System  (HOPS)  and  Error  Subspace  Sta¬ 
tistical  Estimation  (ESSE)  system  were  utilized  in  real-time  to  forecast  physical 
fields  and  uncertainties,  assimilate  various  ocean  measurements  (CTD,  AUVs,  glid¬ 
ers  and  SST  data),  provide  suggestions  for  adaptive  sampling,  and  guide  dynamical 
investigations.  The  qualitative  evaluations  of  the  forecasts  showed  that  many  of  the 
surface  ocean  features  were  predicted,  but  that  their  detailed  positions  and  shapes 
were  less  accurate.  The  root-mean-square  errors  of  the  real-time  forecasts  showed 
that  the  forecasts  had  skill  out  to  2  days.  Mean  one  day  forecast  temperature  RMS 
error  was  0.26°C  less  then  persistence  RMS  error.  Mean  two  day  forecast  tempera¬ 
ture  RMS  error  was  0.13°C  less  then  persistence  RMS  error.  Mean  one  or  two  day 
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salinity  RMS  error  was  0.036PSU  less  than  persistence  RMS  error.  The  real-time 
skill  in  the  surface  was  found  to  be  greater  than  the  skill  at  depth.  Pattern  Correla¬ 
tion  Coefficient  comparisons  showed,  on  average,  greater  skill  than  the  RMS  errors. 
For  simulations  lasting  10  or  more  days,  uncertainties  in  the  boundaries  could  lead 
to  errors  in  the  Monterey  Bay  region. 

Following  the  real-time  experiment,  a  reanalysis  was  performed  in  which  improve¬ 
ments  were  made  in  the  selection  of  model  parameters  and  in  the  open  boundary 
conditions.  The  result  of  the  reanalysis  was  improved  long-term  stability  of  the  sim¬ 
ulations  and  improved  quantitative  skill,  especially  the  skill  in  the  main  thermocline 
(RMS  simulation  error  1°C  less  then  persistence  RMS  error  out  to  5  days).  This 
allowed  for  an  improved  description  of  the  ocean  features.  During  the  experiment 
there  were  two  week-to- 10-day  long  upwelling  events.  Two  types  of  upwelling  events 
were  observed:  one  with  plumes  extending  westward  at  Point  Aho  Nuevo  and  Point 
Sur;  the  other  with  a  thinner  band  of  upwelled  water  parallel  to  the  coast  and 
across  Monterey  Bay.  During  strong  upwelling  events  the  flows  in  the  upper  10-20m 
had  scales  similar  to  atmospheric  scales.  During  relaxation,  kinetic  energy  becomes 
available  and  leads  to  the  development  of  mesoscale  features.  At  100-300m  depths, 
broad  northward  flows  were  observed,  sometimes  with  a  coastal  branch  following 
topographic  features.  An  anticyclone  was  often  observed  in  the  subsurface  fields  in 
the  mouth  of  Monterey  Bay. 


Key  words:  Data  Assimilation,  Oceanography,  Prediction,  Skill,  Uncertainty, 
Upwelling 


1  Introduction 


A  large  Office  of  Naval  Research  (ONR)-sponsored,  multi-institution  coastal 
predictive  skill  exercise,  the  Autonomous  Ocean  Sampling  Network-II  (AOSN- 
II) 1 ,  occurred  in  August  2003  in  the  Monterey  Bay  region  off  central  Califor¬ 
nia.  The  goal  of  this  exercise  was  to  initiate  at-sea  research  of  an  adaptive  ob¬ 
serving  and  prediction  system,  with  the  intent  to  assimilate  various  data  types, 
adapt  the  deployment  of  platforms  and  allow  the  relocation  of  the  system  to 
other  regions  (Bellingham  and  Zhang,  2005).  As  described  in  this  special  issue 
(Ramp  et  ah,  2008;  Curtin  and  Bellingham,  2008;  Davis  et  ah,  2008),  remote 
and  in  situ  sensors  and  platforms  were  employed,  including  gliders  (Davis  et 
ah,  2002;  Rudnick  et  ah,  2004;  Shulman  et  ah,  2008),  drifters  (Chavez  et  ah, 
1997),  moorings,  Autonomous  Underwater  Vehicles  (AUVs),  research  vessels 
and  satellites  and  coastal  radar  (Shulman  and  Paduan,  2008)  for  remote  sens- 

*  Corresponding  author,  phaley@mit.edu 

1  http://www.mbari.org/aosn.  Also  see  Curtin  et  ah,  1989  and  Curtin  and 
Bellingham  2001. 
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ing.  The  collection  of  such  a  data  set,  with  the  intensive  successful  operations 
of  two  fleets  of  gliders  and  low  flying  aircraft  in  coordination  with  vehicle  con¬ 
trol  algorithms  (Leonard  and  Graver,  2001;  Fiorclli  et  ah,  2006;  Leonard  et 
al.,  2007;  Davis  et  al.,  2008)  and  two  modeling  systems,  is  unprecedented  in 
oceanography  (Bellingham,  2006).  The  present  manuscript  reports  the  real¬ 
time  modeling  research  carried  out  with  the  Harvard  Ocean  Prediction  System 
(HOPS;  Robinson  et  al.,  2002).  For  details  of  the  use  of  the  second  modeling 
system  (Rutgers  Ocean  Modeling  System  (ROMS)-West;  Wang  et  al.,  2005; 
Li  et  al.,  2006)  during  AOSN-II,  the  reader  is  referred  to  Chao  et  al.  (2008). 

HOPS  is  an  interdisciplinary  ocean  modeling  system  primarily  designed  for  re¬ 
gional  applications  with  modular  schemes  for  rapid  set-up,  data  assimilation 
and  dynamical  studies  (Robinson  et  al.,  1996;  Robinson,  1999;  Haley  et  al., 
1999).  The  heart  of  HOPS  is  a  physical  primitive-equation  (PE)  dynamical 
model  (Cushman- Roisin,  1994;  Bryan,  1969),  which  is  supported  by  topogra¬ 
phy  conditioning  software,  data  processing  and  gridding  routines,  initialization 
and  assimilation  schemes  (Robinson  et  al.,  1998),  dynamical  studies  schemes, 
and  visualization  software.  During  AOSN-II,  HOPS  assimilated  the  various 
data  types  mentioned  above  and  carried  out  predictions  of  ocean  tempera¬ 
ture,  salinity  and  velocity  fields  over  2-3  days.  ESSE  (Lermusiaux  et  al.,  2002) 
is  a  system  for  the  distributed  prediction  of  oceanic  uncertainties  and  assim¬ 
ilation  of  various  data  types.  ESSE  schemes  (Lermusiaux,  2006)  are  based 
on  an  ensemble  approach  for  the  error  prediction,  on  a  reduction  of  the  er¬ 
ror  space  to  a  dominant  subspace,  and  on  Kalman  Filtering  and  Smoothing 
updates  for  the  assimilation.  During  AOSN-II,  they  were  used  to  carry  out 
ensembles  of  nonlinear  stochastic  forecasts  of  physical  fields  and  uncertainties 
(Lermusiaux,  2006;  Lermusiaux  et  al.,  2006;  Lermusiaux,  2007)  and  to  assim¬ 
ilate  various  data  types.  In  total,  24  sets  of  real-time  nowcasts  and  forecasts 
were  released  from  4  August  to  3  September.  The  forecasts  were  forced  by 
3km  and  hourly  ocean-atmosphere  flux  predictions  obtained  from  the  U.S. 
Navy’s  operational  Coupled  Ocean/ Atmosphere  Mesoscale  Prediction  System 
(COAMPS,  Doyle  et  al.,  2008).  The  ocean  observations  collected  were  qual¬ 
ity  controlled  and  mapped  by  objective  analysis  daily.  Boundary  conditions 
and  model  parameters  for  atmospheric  forcing  were  calibrated  and  modified 
in  real-time  to  adapt  to  the  evolving  ocean  conditions.  Products  were  dis¬ 
tributed  on  the  Web,  oceanic  features  and  uncertainty  fields  were  described 
and  adaptive  sampling  plans  were  provided  on  a  daily  basis. 

The  main  features  of  the  AOSN-II  region  are  diagrammed  in  Fig.  1,  along 
with  the  chosen  HOPS  modeling  domains.  These  features  were  obtained  from 
a  study  of  the  literature  and  from  the  results  of  the  AOSN-II  exercise.  They 
are  the:  (1)  Upwelling  centers  at  Point  Aho  Nuevo  (AN)  and  Point  Sur  (PS), 
with  the  upwelled  water  advected  equatorward  and  seaward  (Huyer,  1983; 
Rosenfeld  et  al.,  1994;  Traganza  et  al.,  1981;  Ramp  et  al.,  2005);  (2)  Coastal 
current,  eddies,  squirts  and  filaments  including  the  upwelling- induced  jets  and 
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high  (sub)-mesoscale  variability  in  the  Coastal  Transition  Zone  (Brink  et  ah, 
1991;  Huyer  et  ah,  1991;  Kosro  et  ah,  2002;  Strub  et  ah,  1991;  Ramp  et  ah, 
1991;  Hayward  and  Mantyla,  1990;  Chavez  et  ah,  1991;  Chavez  and  Collins, 
2000);  (3)  California  Undercurrent  (CUC),  which  is  a  poleward  flow/jet  often 
found  10-100km  offshore  and  at  50-300m  depths  (Collins  et  ah,  1996;  Garfield 
et  ah,  1999;  Pierce  et  ah,  2000;  Wickham,  1975);  (4)  California  Current  (CC), 
a  broad  southward  flow,  often  found  100-1350km  offshore  and  at  0-500m  depth 
(Collins  et  ah,  2003;  Huyer  et  ah,  1998;  Hickey,  1998;  Marchesiello  et  ah,  2003; 
Lynn  and  Simpson,  1987;  Strub  and  James,  2000). 

In  the  Monterey  Bay  region,  the  major  external  drivers  of  the  surface  ocean 
dynamics  are  atmospheric  fluxes,  especially  the  wind  stress.  In  fact,  these 
forcings  define  two  different  ocean  states;  classically  known  as  upwelling  and 
relaxation  states  (Rosenfeld  et  ah,  1994).  To  provide  a  time  line  of  these  dif¬ 
ferent  events  during  AOSN-II,  Figure  2  shows  the  measured  winds  at  the  Ml 
and  M2  moorings.  In  the  period  of  August  8-17,  both  moorings  show  a  sus¬ 
tained  upwelling  event,  with  a  brief  weakening  around  August  15.  Between 
August  17-21,  the  more  offshore  M2  shows  a  continuation  of  upwelling  favor¬ 
able  winds  while  Ml  indicates  weakening  and  less  uniformity  in  direction.  The 
two  moorings  indicate  a  relaxation  event  between  August  21-24.  After,  during 
August  24-27,  M2  shows  the  start  of  another  upwelling  event  but  Ml,  closer 
to  the  shore,  again  shows  less  uniformity  in  direction.  Between  August  27  and 
September  1,  the  winds  are  generally  in  an  upwelling  favorable  direction,  al¬ 
though  the  strength  is  falling  off  after  August  29.  Both  moorings  indicate  a 
relaxation  event  between  September  2-5. 

Figure  3  shows  the  positions  of  the  main  sources  of  in  situ  data  used  in  this 
study.  Three  comprehensive  CTD  surveys  were  made  by  the  R/V  Pt.  Sur  (2- 
6  Aug;  21-25  Aug  and  3-6  Sep  2003).  Each  survey  consisted  of  58-69  CTD 
stations  (69  for  the  first  2,  58  for  the  third)  and,  for  the  purposes  of  this 
study,  each  defined  a  natural  time-window  for  constructing  synoptic  estimates 
of  the  ocean  state.  Also  shown  are  the  positions  of  the  pseudo-profiles  from 
the  gliders.  Scripps  (SIO)  maintained  5  gliders  (Sherman  et  ah,  2001)  in  the 
period  23  Jul  -  2  Sep,  2003;  traveling  yo-yo  paths  between  the  surface  and 
400m.  Woods  Hole  (WHOI)  maintained  10  gliders  in  the  period  21  Jul  -  31 
Aug,  2003,  sampling  between  the  surface  and  200m.  These  gliders  provided 
the  bulk  of  the  in  situ  ocean  data  (for  more  details  of  the  AOSN-II  glider 
operations,  see  section  4  of  Rudnick  et  ah,  2004). 

The  two  domains  represented  in  Fig.  1  were  designed  to  support  2-way  nested 
simulations.  The  innermost  domain,  referred  to  as  the  “Data  Domain”,  con¬ 
tains  the  region  of  primary  sampling  around  Monterey  Bay,  as  well  as  the 
upwelling  centers  at  Point  Ano  Nuevo  and  Point  Sur.  The  larger  domain,  re¬ 
ferred  to  as  the  “Off  Shore  Domain”,  provides  a  buffer  for  the  Data  Domain 
while  capturing  the  longest  expected  plumes.  The  term  “Data  domain”  was 
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chosen  based  on  the  AOSN-II  sampling  plans  (Fig  3). 


During  the  exercise,  the  HOPS  modeling  was  carried  out  in  three  different 
configurations.  The  first  was  the  “Stand-alone  Data  Domain”.  In  this  case 
the  PE  model  was  conhgured  to  run  in  the  Data  Domain  with  open  boundary 
condition  schemes  (§2.3.4)  utilized  at  the  external  boundaries.  Data  assimila¬ 
tion  was  done  using  the  Optimal  Interpolation  (01)  method  (Lozano  et  al., 
1996;  Lermusiaux,  1999).  This  configuration  is  always  tried  in  a  new  region,  so 
as  to  provide  a  modeling  benchmark.  Its  relative  simplicity  facilitates  a  rapid 
search  for  a  stable  regime.  The  second  configuration  also  used  01  assimilation 
but  was  set  for  a  2-way  nested  communication  (Sloan,  1996;  Spall  and  Hol¬ 
land,  1991;  Fox  and  Maskcll,  1995)  between  one  model  code  set  up  for  the 
Data  Domain  and  another  for  the  Off  Shore  Domain  (see  App.  I  for  details). 
Theoretically,  this  provides  a  superior  representation  of  the  larger-scale  lateral 
forcing  on  the  smaller  domain  while  simultaneously  reducing  boundary  noise. 
However,  if  these  larger-scale  estimates  are  not  sufficiently  accurate,  nesting 
can  damage  the  fields  in  the  smaller  data-driven  domain.  During  AOSN-II, 
there  were  no  significant  in  situ  data  outside  of  the  Data  Domain  and  the 
larger-scale  synoptic  estimates  were  thus  of  limited  reliability.  It  was  only 
after  shifting  the  GDEM 2  climatology  (Davis  et  ah,  1986)  to  the  observed 
AOSN-II  mean  that  the  larger-scale  Off  Shore  Domain  was  of  some  use.  The 
third  configuration  was  the  “ESSE”  configuration.  This  was  essentially  the 
“Stand-alone  Data  Domain”  with  assimilation  being  done  via  the  ESSE.  This 
configuration  had  a  superior  data  assimilation  methodology  and  also  produced 
estimates  of  the  forecast  error.  The  forecasted  error  fields  are  not  described 
in  this  manuscript.  Although  3  model  configurations  were  used,  on  a  given 
day  the  issued  forecast  products  came  from  a  single  configuration  (except  for 
Aug  12-13).  One  reason  for  this  was  that  it  was  regarded  by  the  other  mem¬ 
bers  of  the  AOSN-II  team  as  confusing  on  the  two  days  for  which  we  issued 
from  multiple  sources.  More  importantly,  we  did  not  always  have  simultane¬ 
ously  useful  forecasts  from  all  configurations.  As  we  tried  new  initialization 
schemes  (§2.3.2)  we  would  first  find  stable  parameter  regimes  in  the  stand¬ 
alone  01  configuration  before  attempting  a  nested  configuration.  At  one  point, 
the  stand-alone  ESSE  configuration  was  being  used  for  extensive  parameter 
tuning,  including  new  boundary  conditions  (App.  II),  following  which  only 
ESSE  forecasts  were  available  while  the  other  configurations  caught  up. 

In  what  follows,  the  real-time  approach  and  modeling  system  components 
are  first  summarized  (Sect  2).  The  real-time  results  are  presented  in  Sect  3, 
including  a  description  of  the  general  circulation  and  hydrodynamics  estimates 
and  evaluations,  both  qualitative  and  quantitative,  of  the  forecast  skills  and 
features.  The  re-analysis  methods  and  results  are  described  in  Sect.  4  and  5, 
respectively.  Results  include  a  synthesis  of  ocean  features  and  properties  for 


2  https : // 128 . 160 . 23 . 42/gdemv/gdemv . html 
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the  region  and  a  evaluation  of  the  re- analysis  skill.  Conclusions  are  in  Sect.  6. 


2  Real-Time  Approach,  Constraints  and  Modeling  System  Com¬ 
ponents 


2. 1  Approach 


The  chosen  scientific  approach  sets  the  forecast  methodologies  and  operational 
procedures.  Our  real-time  approach  is  to:  (i)  initialize  from  a  synoptic  survey 
so  as  to  capture  the  background  ocean  fields  from  the  start;  (ii)  subsequently 
utilize  as  much  ocean  data  as  possible,  both  for  skill  evaluation  and  assimi¬ 
lation,  (iii)  keep  some  data  out  of  the  assimilation  for  sustained  independent 
evaluations  and,  (iv)  allow  human  intervention  and  carry  out  real-time  sci¬ 
entific  studies  of  the  outputs  prior  to  product  dissemination.  Detailed  exam¬ 
ination  of  each  data-driven  forecast  output  is  not  a  guarantee  of  successful 
prediction  but  it  is  nonetheless  essential.  It  involves  studying,  at  a  minimum, 
the  evolution  of  the  ocean  fields  for  each  forecast  candidate  at  the  depths 
and  along  the  vertical  sections  that  are  dynamically  key.  For  AOSN-II,  this 
was  mainly  human-based  but  research  is  underway  for  autonomous  scientific 
evaluations  (Lermusiaux,  2007). 


2.2  Constraints  and  Operations 


A  main  constraint  on  the  above  choices  and  approach  is  the  ability  to  pro¬ 
vide  daily  forecast  products  in  a  timely  manner.  This  “timeliness”  constraint 
is  set  by  the  available  personnel  and  computing  power  3 .  For  AOSN-II,  the 
chosen  procedure  was  to  prepare  a  set  of  candidate  forecasts  each  morning 
so  as  to  include  the  most  recent  data  and  atmospheric  forcings.  The  main 
timeliness  constraint  was  that  the  forecasts  needed  to  be  finished,  examined 
and  have  their  products  generated  by  late  afternoon  for  the  AOSN-II  team 
teleconferences. 

The  canonical  day’s  operations  then  proceeded  as  follows  (see  also  fig.  4).  First, 
the  most  recent  data  and  atmospheric  forcings  were  gathered  and  processed  for 
use  in  the  HOPS  PE  model  and  ESSE  ensemble  predictions.  An  initialization 
was  generated  from  the  previous  day’s  simulation,  from  the  fields  correspond¬ 
ing  to  0000Z  on  the  previous  day.  To  calibrate  numerical  and  dynamical  model 

3  For  AOSN-II,  we  utilized  a  set  of  10  to  18  distributed  processors  4  to  6  years  old 
and  with  200  to  400Mhz  clock  speed 


6 


parameters  to  the  region  and  data,  several  PE  simulations  would  then  be  run, 
using  the  updated  atmospheric  forcings  and  assimilating  the  new  data  but 
with  different  model  parameters.  These  simulations  were  then  visually  com¬ 
pared  to  available  data  (CTD,  glider,  aircraft  and  satellite  SST,  CODAR) 
assessing  the  evolution  of  the  main  features  (depths  of  the  thermocline  and 
surface  mixed  layer,  extent  of  any  upwclling,  and  general  circulation  in  Mon¬ 
terey  Bay).  The  simulations  were  also  inspected  for  the  presence  and  extent 
of  numerical  error  (e.g.  errors  originating  at  the  open  boundaries).  The  simu¬ 
lation  with  the  best  data  match  and  minimal  numerical  error  would  then  be 
selected  and  the  corresponding  forecast  products  generated  and  disseminated. 
There  were  2  paths  for  the  dissemination  of  HOPS  products.  The  Erst  was 
our  local  web  page  4  ,  the  second  was  the  MBARI  server  5  .  Adaptive  sampling 
recommendations  were  also  provided  to  the  AOSN-II  team  on  a  daily  basis 
(see  RTOC  links  in  local  web  site).  Their  heuristic  objectives  were  (i)  to  sam¬ 
ple  interesting  features  of  the  predicted  dynamics,  (ii)  to  reduce  the  predicted 
uncertainties  of  the  ocean  held  estimates  or  (iii)  to  maintain  coverage  over  the 
whole  ocean  region.  Results  from  this  adaptive  sampling  have  been  previously 
published  by  Lermusiaux  (2007). 


2.3  Modeling  System  Components 


In  real-time  regional  ocean  predictions,  the  choices  made  in  setting  up  the 
modeling  system  for  the  dynamics  of  interest  are  essential.  This  is  because  a 
modeling  system  consists  of  a  large  number  of  components  that  need  to  be  ad¬ 
equately  selected  given  the  constraints  of  the  specific  prediction  exercise.  The 
components  of  modeling  systems  include:  the  bathymetry,  the  domains  (nest¬ 
ing,  stand-alone,  etc),  the  numerical  properties  (scheme,  grids,  vertical  and 
horizontal  resolution,  etc),  the  data  utilized,  the  initial  conditions,  the  open 
and  land  boundary  conditions,  the  forcings  including  atmospheric  fluxes,  tides 
and  rivers,  the  data  assimilation  scheme,  the  chosen  model  dynamics,  the  dy¬ 
namical  parameterizations  and  finally,  the  dynamical  and  computational  pa¬ 
rameters.  Some  of  these  components  are  summarized  next.  Additional  details 
are  given  in  appendix  I. 


2.3.1  Data 

The  data  utilized  for  evaluation  and  assimilation  consisted  of  measurements 
collected  by  gliders,  aircraft,  AUVs  and  ships  (Ramp  et  al.,  2008).  The  syn¬ 
optic  oceanographic  data  were  evaluated,  processed  and  objectively  analyzed 
on  a  daily  basis.  The  data  were  first  visually  inspected  on  an  instrument  by 

1  http :  //people .  deas  .harvard.  edu/^leslie/AOSNII/ 

5  http :  //www.mbari  .  org/aosn/A0SN_MB2003table . htm 
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instrument  basis  to  ensure  internal  consistency.  In  some  cases,  filtering  was 
applied.  After  this  first  pass  of  quality  control,  data  from  different  instruments 
were  visually  compared  to  each  other.  Outliers  were  identified  and  removed, 
retaining  as  much  data  as  possible.  Once  the  data  had  passed  this  quality  con¬ 
trol  procedure,  they  were  objectively  analyzed  (Bretherton  et  ah,  1976;  Carter 
and  Robinson,  1987). 


2.3.2  Initialization 

For  the  AOSN-II  region,  the  initial  conditions  were  found  to  set  the  quality 
of  the  ocean  predictions  over  days  to  weeks.  This  is,  in  a  large  part,  related 
to  the  time  required  to  modify  the  mean  background  state  below  the  surface 
mixing  layer  and  over  the  larger  domain  shown  on  Fig.  1.  During  the  AOSN- 
II  exercise,  the  initialization  fields  we  utilized  changed  with  time,  according 
to  the  three  R/V  Pt.  Sur  surveys.  Prior  to  August  7,  the  first  R/V  Pt.  Sur 
survey  data  were  not  yet  completed  and  processed.  The  forecasts  issued  in 
those  times  were  initialized  with  July  2003  data  combined  with  historical 
synoptic  surveys.  Three  other  initial  conditions  based  on  climatological  data, 
El  Nino  data  and  La  Nina  data  were  also  utilized  but  none  matched  the 
synoptic  data  as  it  was  collected.  The  July  2003  data  set  was  a  better  fit,  but 
still  not  a  very  good  one.  In  the  period  August  7-26,  the  processed  first  survey 
was  available.  The  forecasts  were  then  either  directly  re-initialized  with  this 
first  survey  or  were  restarts  from  forecasts  that,  ultimately,  traced  back  to  this 
initialization  survey.  After  August  26,  the  second  survey  was  processed  and  the 
forecasts  were  again  either  initialized  from  this  second  survey  or  from  restarts 
that  traced  back  to  it.  This  applies  to  both  the  fields  and  to  the  uncertainty 
computations  using  ESSE.  This  re- initialization  was  useful  to  reduce  the  effects 
of  uncertainties  in  open  boundary  conditions  which  had  been  compounded  for 
3  weeks  of  forecasting  from  the  previous  survey. 


2.3.3  Forcings 

The  atmospheric  forcings  for  our  AOSN-II  HOPS  simulations  were  generated 
based  on  the  Navy’s  operational  Coupled  Ocean/ Atmosphere  Mesoscale  Pre¬ 
diction  System  (COAMPS)  (Hodur,  1997)  analysis  and  72  hour  forecast  fields, 
received  on  a  twice  daily  basis  (Doyle  et  ah,  2008).  The  COAMPS  fields  were 
available  on  four  different  resolution  lambert  conformal  projection  grids.  The 
nominal  0.03°  resolution  fields  were  used  to  generate  forcings  for  the  Data  Do¬ 
main.  The  nominal  0.1°  fields  were  used  to  generate  forcings  for  the  Offshore 
domain.  Overall,  it  was  found  that  these  forcings  were  of  a  quality  superior 
to  anything  we  had  utilized  before.  We  also  found  that  the  0.03°  resolution 
forcing  was  necessary  for  ocean  predictions  of  upwelling  events  at  AN  and  PS. 
Finally,  our  AOSN-II  simulations  were  not  forced  by  tides  nor  by  river  inputs 


(Rosenfeld  et  al.,  2008;  Wang  et  al.,  2008). 


2.3.4  Parameterizations  and  Parameters 

In  our  general  approach  (§2.1),  we  usually  test  our  system  with  relevant  his¬ 
torical  synoptic  datasets  before  the  real-time  operation  so  as  to  estimate  a 
range  of  adequate  parameter  values.  In  the  six  months  prior  to  AOSN-II, 
we  were  provided  with  several  datasets:  June- August  2000  ICON  primitive- 
equation  model  fields  and  atmospheric  fluxes  from  Igor  Shulman  (Shulman 
et  ah,  2002);  a  May-June  1989  synoptic  survey  from  Leslie  Rosenfeld;  as  well 
as  NODC  and  CalCOFI  (Scripps,  1999)  data.  Using  these  data  with  synop¬ 
tic  winds  and  process-oriented  idealized  winds  (upwelling- relaxation  cycles), 
we  chose  an  initial  set  of  dynamical  and  computational  parameters  for  both 
numerical  stability  and  dynamical  response. 

The  parameterizations  and  parameter  value  ranges  that  were  utilized  dur¬ 
ing  AOSN-II  were  as  follows.  All  simulations  were  made  with  a  300  second 
time  step.  The  lateral  open  boundary  forcing  was  parameterized  with  the  im¬ 
plicit  Orlanski  radiation  condition  (Orlanski ,  1976)  for  the  tracer,  velocity  and 
transport  streamfunction  variables.  For  the  rate  of  change  of  barotropic  vortic- 
ity,  the  boundary  condition  most  often  used  was  a  CFvN  condition  (Charney 
et  ah,  1950),  recast  for  the  PE  by  Spall  and  Robinson  (1989).  In  the  last  third 
of  the  experiment,  the  Orlanski  condition  was  used  on  vorticity.  The  hori¬ 
zontal  sub-gridscale  was  parameterized  by  a  Shapiro  filter  (Shapiro,  1970):  a 
fourth  order  filter  was  used  for  tracers  and  momentum,  and  a  (stronger)  sec¬ 
ond  order  filter  for  the  rate  of  change  of  barotropic  vorticity.  In  the  vertical, 
the  background  sub-gridscale  was  parameterized  by  a  second  order  diffusion 
term  where  the  coefficients  were  functions  of  the  Richardson  number,  as  in 
(Pacanowski  and  Philander,  1981).  The  peak  viscosity/ diffusivity  was  set  to 
50cm2/s.  The  background  viscosity  and  diffusivity  were  set  in  the  range  0.1- 
0.2cm2/s  and  0.01-0.02cm2/s,  respectively.  Near  the  surface,  vertical  mixing 
(see  Lermusiaux,  2001)  is  increased  to  at  least  30cm2/s  for  viscosity  and  5- 

6cm2 /s  for  diffusivity,  so  as  to  reflect  wind  forcing.  The  depth  to  which  these 

1 

mixing  bounds  is  applied  is  proportional  to  (||r||  /po)2  //  where  r  is  the  wind 
stress,  po  is  the  mean  density  of  seawater  and  /  is  the  Coriolis  factor.  The 
nondimensional  constant  of  proportionality  was  set  in  the  range  0.15-0.22. 
This  wind  mixing  depth  was  never  allowed  to  be  shallower  than  lrn  or  ex¬ 
ceed  40m.  At  depths  where  the  water  column  is  gravitationally  unstable,  the 
vertical  viscosity/ diffusivity  were  bounded  from  below  by  50cm2/s.  Finally,  a 
Rayleigh  friction  parameterization  of  drag  induced  by  the  bottom  or  coasts 
(Lermusiaux,  1997)  was  used.  In  both  cases,  the  temporal  decay  scale  was 
3600s  and  the  spatial  decay  was  1.5  (model  levels  from  the  bottom  or  grid 
points  from  the  nearest  coast,  respectively). 


9 


2.3.5  Data  Assimilation 


The  assimilation  procedure  was  based  on  a  daily  cycle 6 .  Once  each  simu¬ 
lation  day,  data  were  assimilated  at  a  nominal  assimilation  time  of  1200Z. 
Temperature  and  salinity  data  within  ±18hr  of  that  time  were  gathered  and 
objectively  analyzed  with  the  assimilation  scales  in  table  1.  These  fields  were 
then  processed  for  the  terrain-following  coordinates,  with  a  geostrophic  in¬ 
ternal  mode  velocity.  These  temperature,  salinity  and  internal  velocity  fields 
were  then  assimilated  into  the  model,  according  to  their  associated  error  fields. 
Rather  than  “shocking”  the  system  by  simply  inserting  the  assimilation  fields 
at  the  desired  times,  the  fields  are  first  assimilated  earlier  with  down-weighted 
assimilation  weights  (“ramping  in”).  The  schedule  used,  along  with  the  down¬ 
weighting  weights,  are  given  in  table  2.  As  table  2  indicates,  we  found  out  that 
the  weights  for  assimilating  the  internal  mode  geostrophic  velocity  needed  to 
be  smaller  than  those  for  temperature  and  salinity,  so  as  to  not  locally  over¬ 
write  the  PE  physics  with  geostrophy.  A  weak  assimilation  after  the  center 
time  (a  “ramp  out”  cycle)  was  also  included.  Note  that  prior  to  August  19, 
the  “ramp  out”  cycle  was  not  used.  Based  on  forecast  evaluations,  the  weights 
for  assimilating  internal  mode  velocity  were  further  decreased  by  a  factor  of 
2,  on  August  19-21,  and  28-30. 


3  Real-Time  Results 

3.1  General  Circulation  and  Hydrodynamics  Estimates 


From  August  4  -  September  3,  2003,  forecasts  were  issued  on  23  days.  Figure  5 
shows  the  temperature  and  velocity  fields  for  8  Aug  2003,  which  are  illustrative 
of  the  state  of  the  general  circulation  and  hydrographic  features  during  the 
entire  period.  This  state  was  found  to  be  in  stronger-than-usual  upwclling 
favorable  conditions  (Ramp  et  ah,  2008).  Starting  with  the  30m  fields,  the 
forecasts  show  a  general  cyclonic  circulation  in  Monterey  Bay,  a  situation  we 
have  found  to  occur  during  upwelling  conditions.  For  significant  upwclling 
conditions,  a  southward  coastal  current  was  often  observed  flowing  along  the 
shelf  and  across  the  mouth  of  Monterey  Bay.  Offshore,  we  observed  an  inflow 
in  the  western  boundary,  bending  north  and  providing  a  general  northward 
flow.  This  is  indicative  of  a  surfacing  undercurrent  also  referred  to  as  the 
Davidson  current  (Brink  et  ah,  1991).  Inshore  of  this  flow,  but  off  the  shelf, 
an  anticyclonic  eddy  was  usually  present.  The  exact  size,  shape  and  position 

6  A  finding  from  AOSN-II  is  that  the  assimilation  should  account  for  faster  diurnal 
and  tidal  cycles.  For  subsequent  experiments,  twice  daily  assimilation  cycles  were 
used. 
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of  this  eddy  was  found  to  be  fairly  variable.  The  situation  at  the  surface  is 
very  similar  except  during  the  periods  of  strong  upwelling  favorable  winds. 
In  these  periods,  winds  dominate,  resulting  in  broad  generally  southward  flow 
across  the  whole  domain.  The  surface  also  shows  a  band  of  warmer  water 
lining  the  coast  of  Monterey  Bay,  occasionally  protruding  out  the  northern 
edge  of  the  mouth  of  the  bay  (not  on  Aug  8).  This  feature  was  sometimes  in 
SST  images.  Our  forecasts  at  times  over-estimated  the  extent  of  the  northward 
protrusion,  in  part  due  to  fog  conditions  perhaps  not  fully  represented  in  the 
atmospheric  forcing  fluxes  but  also  in  part  due  to  flow  fields  that  were  too 
geostrophic.  At  200m,  there  is  a  general  northward  flow  originating  from  the 
western  boundary,  sometimes  accompanied  by  a  northward  branch  along  the 
slope,  characteristic  of  the  undercurrent.  The  eddy  opposite  Monterey  Bay  is 
still  visible  at  this  depth,  sometimes  even  more  so  than  at  shallower  depths. 


Once  the  larger-scale  initial  conditions  were  estimated  by  the  first  R/V  Pt. 
Sur  survey,  several  uncertainties  still  remained  in  our  (sub)-mesoscale  predic¬ 
tions.  This  can  be  illustrated  by  predictions  made  by  the  three  different  PE 
configurations.  Consider  the  forecasts  that  were  issued  on  August  13th  (fig¬ 
ure  6).  An  important  difference  can  be  seen  at  200m  (bottom  row).  Although 
the  three  cases  show  general  northward  flow,  the  stand-alone  case  has  2  dis¬ 
tinct  branches.  The  first  enters  from  the  western  boundary  and  turns  north. 
The  second  enters  from  the  southern  boundary  by  Point  Sur  and  flows  north 
along  the  slope.  This  second  branch  pushes  an  eddy,  originally  opposite  Mon¬ 
terey  Bay,  northward  along  the  slope  until,  by  the  August  13,  it  is  opposite 
point  Ano  Nuevo.  In  contrast,  the  nested  configuration  has  a  large,  anticy- 
clonic  half  eddy  on  the  southern  boundary.  This  boundary  eddy  appears  to 
choke  off  the  along  slope  branch  of  the  northward  flow.  As  a  result,  the  eddy 
that  was  originally  opposite  Monterey  Bay  does  not  advect  nearly  so  far  north 
by  August  13.  The  inflow  from  the  western  boundary  is  broader  and  further 
south.  It  still  bends  north  in  a  broad  northward  flow  that  fills  most  of  the 
domain.  The  ESSE  forecast  is  more  similar  to  the  stand-alone  case,  in  that 
it  has  2  branches.  However,  the  inflow  from  the  west  is  much  broader  for  the 
ESSE  forecast  than  the  stand  alone  01  forecast,  covering  the  northern  2/3 
of  the  western  boundary.  Overall,  the  open  boundary  behaves  better  in  the 
ESSE  forecast  than  in  the  stand  alone  OI  forecast.  Since  the  initial  conditions 
were  similar  in  the  144  km  by  124.5  km  Data  Domain,  a  conclusion  is  that 
over  10  days  of  data-driven  forecasting,  the  open  boundary  forcing  has  a  sig¬ 
nificant  impact  on  the  Monterey  Bay  region.  Uncertainties  also  exist  in  the 
atmospheric  forcing.  Although  the  general  upwclling/relaxation  cycles  were 
well  represented,  the  atmospheric  uncertainties  would  primarily  affect  the  de¬ 
tails  of  the  fields  in  the  upper  layers  of  the  forecasts  (e.g.  the  warm  protrusion 
of  the  preceding  paragraph). 
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3.2  Qualitative  evaluation  of  forecast  circulation  and  hydrographic  features 


As  the  AOSN-II  predictive  skill  exercise  proceeded,  qualitative  assessments 
were  made  with  whatever  data  and  images  were  available  in  real-time  and 
reported  7 .  A  few  of  these  real-time  assessments  are  reproduced  in  figures  7  and 
8,  focusing  on  comparisons  to  SST  and  CODAR  data.  Comparing  the  HOPS 
surface  temperature  1  day  forecast  to  satellite  SST  on  12  Aug  2003  (left  panels 
Fig.  7),  the  model  reproduces  the  general  feature  of  an  upwclling  plume  from 
the  north  being  advected  across  the  mouth  of  Monterey  Bay  and  joining  with 
the  upwelling  off  Point  Sur.  The  model  seems  to  advect  more  of  the  surface 
material  to  the  south  and  fails  to  capture  a  westward  extension  of  the  plume 
by  the  northern  edge  of  Monterey  Bay.  This  is  indicative  of  inaccuracies  in 
the  surface  boundary  layer  parameters  (they  were  in  fact  changed  in  real-time 
that  week:  see  Lermusiaux,  2007)  as  well  as  inaccuracies  in  the  atmospheric 
forcing  fluxes.  Surface  properties  are  very  sensitive  to  such  parameters:  for 
example,  an  error  of  25  percent  over  a  surface  flow  of  60cm/s  leads  to  an  error 
of  13  km  per  day.  On  13  Aug  2003  (middle  panels  of  Fig  7),  a  comparison 
of  the  surface  temperature  to  aircraft  SST  still  shows  the  cold  plume  across 
the  mouth  of  the  bay.  Both  the  model  surface  temperature  and  the  aircraft 
SST  display  warmer  waters  right  at  the  coast  of  Monterey  Bay,  although  the 
model  fields  have  these  waters  warmer  and  have  a  little  warm  water  plume 
extending  to  the  west  out  of  the  northern  edge  of  Monterey  Bay  (lack  of  fog 
in  atmospheric  forcing  and  too  geostrophic  flow  in  the  ocean  model).  Looking 
at  the  unobstructed  portion  of  the  satellite  SST  image  from  21  Aug  2003 
(right  panels  Fig.  7)  shows  an  elongated  plume  extending  westward  out  of 
the  southern  edge  of  Monterey  Bay,  an  upwclling  center  over  Point  Sur  and 
a  wedge  of  warmer  water  separating  them.  The  one  day  forecast  for  21  Aug 
2003  captures  the  two  cold  pools  separated  by  a  warmer  wedge.  However, 
the  westward  extension  of  the  plume  is  too  short  and  the  upwelling  center  off 
Point  Sur  is  displaced  slightly  to  the  south.  In  conclusion,  the  predicted  surface 
temperature  fields  contain  the  dominant  features  but  some  of  the  smaller  scales 
and  more  nonlinear  phenomena  were  not  well  captured,  especially  to  the  west 
and  south  of  the  Bay,  where  observations  were  limited. 

Figure  8  shows  the  real-time  comparisons  between  the  available  CODAR  im¬ 
ages  and  the  HOPS  forecasts.  Note  that  both  instantaneous  and  daily-averaged 
CODAR  images  were  used.  This  is  because  our  primary  concern  was  the  gen¬ 
eral  meso-scale  circulation,  which  could  be  inferred  from  either  type  of  image. 
The  19  Aug  2003  nowcast  surface  velocities  show  good  qualitative  agreement 
with  the  mapped  CODAR  data  (left  panels  of  Fig.  8).  Both  exhibit  a  south¬ 
ward  flow  across  the  mouth  of  Monterey  Bay  and  a  general  cyclonic  circulation 
within  the  bay.  By  24  Aug  (middle  panels  of  Fig.  8),  after  the  relaxation,  the 


7  See  the  RTOC  presentations  in  the  local  web  site  (§2.2) 
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former  cyclonic  Bay  circulation  is  no  longer  evident  in  either  the  CODAR  nor 
the  HOPS  fields.  The  inflow  in  the  Bay  is  captured  well,  as  is  the  quiet  north¬ 
east  corner.  The  anti-cyclonic  feature  outside  of  the  bay  is  not  captured.  By  30 
Aug  (right  panels  of  Fig.  8),  CODAR  again  shows  a  cyclonic  circulation  inside 
Monterey  Bay,  along  with  a  southward  flow  at  the  southern  end  of  Monterey 
Bay  and  a  meandering  southward  flow  offshore  of  Monterey  Bay.  All  are  well 
represented  in  the  1  day  forecast  for  30  Aug.  In  conclusion,  the  assimilation  of 
high-resolution  WHOI  glider  data  from  the  northern  portion  of  the  Bay  (see 
Fig.  3)  leads  to  good  predictions  of  the  surface  currents.  This  was  encouraging 
since  the  ocean  model  utilized  for  AOSN-II  did  not  yet  contain  tidal  effects, 
which  have  non-zero  mean  effects  on  the  hydrography  and  circulation. 

Many  other  evaluations  were  carried  out.  One  interesting  one  was  carried  by  J. 
Marsden  and  his  team  at  the  California  Institute  of  Technology  (Chang,  2003). 
They  used  the  HOPS  forecast  velocity  fields  to  compute  Lagrangian  Coherent 
Structures  (LCS)  (Haller,  2002;  Shadden  et  al.,  2005,  2008;  Lermusiaux  and 
Lekien,  2005;  Lermusiaux  et  al.,  2006;  Lermusiaux,  2006)  in  real-time  (using 
ManGen  8  )  and  used  these  structures  to  predict  the  track  of  a  drifter  launched 
in  Monterey  Bay  by  F.  Chavez  at  MBARI.  The  drifter  always  stayed  on  the 
same  side  of  the  predicted  LCS  separatrix,  parallel  to  the  mouth  of  the  Bay  at 
Point  Pinos.  This  demonstrated  that  LCS  are  useful  to  characterize  hyperbolic 
ocean  features  and  that  the  HOPS  surface  current  predictions  were  accurate 
enough  to  estimate  the  LCS.  Admittedly,  this  was  only  a  single  event  and 
therefore  of  limited  probative  value.  However  it  is  usually  the  case  that  real¬ 
time  operations  have  only  limited  (or  even  singular)  opportunities  for  ancillary 
tests  of  the  forecasting  system. 


3.3  Quantitative  evaluation  of  forecast  fields 


To  assess  the  predictive  skill  in  the  real-time  forecasts,  two  skill  metrics  are 
introduced.  The  first  is  the  root-mean-square  error.  Denoting  the  set  of  fore¬ 
cast  values  Tf  and  the  corresponding  set  of  observed  values  T0,  then  the 
root-mean-square  forecast  error,  TrmsForecast ,  is  simply  the  root  mean  square 
of  the  difference  between  Tf  and  T0.  Similarly  the  RMS  persistence  error, 
Trmspersistence,  is  obtained  by  replacing  the  forecast  values  with  the  initializa¬ 
tion  values  T;  in  the  root-mean-square  difference.  A  forecast  is  said  to  have 
skill  when  the  RMS  forecast  error  is  less  than  the  RMS  persistence  error.  A 
large  number  of  skill  studies  were  carried  out.  The  illustrations  shown  below 
summarize  the  results. 

Figures  9-11  show  differences  between  the  RMS  persistence  errors  and  the 
RMS  forecast  errors  for  the  one  and  two  day  forecasts  of  temperature  and 

8  http://www.mangen.info/  see  also  Lekien  et  al.  (2005);  Shadden  et  al.  (2007) 
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salinity  issued  in  real-time.  Each  forecast  is  compared  to  data  not  assimilated 
into  the  forecast.  This  data  is  in  the  form  of  assimilation  fields  that  were 
made  after  the  forecast  was  issued,  using  data  not  available  to  the  forecast. 
The  comparisons  were  restricted  to  the  regions  where  the  nondimensional 
observation  OA  error  is  less  than  or  equal  to  0.25.  The  RMSEs  are  computed 
at  2m  (surface),  10m  (base  mixed  layer/  top  of  thermocline),  30m  (middle  of 
thermocline),  150m  (California  Under  Current),  300m  (deep)  and  the  mean 
of  all  the  selected  points  on  those  levels.  For  a  given  forecast,  persistence  is 
defined  as  the  particular  initialization  (see  Sect  2.3)  from  which  that  forecast 
descends.  That  means  that  the  comparison  is  to  one  of  three  sets  of  fields, 
depending  on  whether  the  simulation  was  made  prior  to  Aug  7,  between  Aug 
7-26  or  after  Aug  26. 

Figure  9  displays  the  temperature  (in  C)  and  salinity  (in  PSU)  differences 
for  the  means  at  1  and  2  day  forecasts.  On  average,  both  temperature  and 
salinity  forecasts  have  skill  out  to  2  days.  The  average  1  day  temperature 
forecast  shows  greater  skill  than  the  2  day  forecast  (0.26  mean  difference  vs 
0.13),  whereas  the  salinity  1  and  2  day  forecasts  show  equal  skill  (0.036  mean 
difference).  The  1  day  forecast  mean  temperature  differences  show  a  correla¬ 
tion  with  the  wind,  having  skill  during  the  upwelling  period  of  August  8-20 
and  a  loss  of  skill  during  the  relaxation  period  of  August  21-24.  A  important 
contributing  factor  to  this  is  the  fact  that  the  persistence  being  used  in  both 
periods  comes  from  the  data  during  the  first  R/V  Pt.  Sur  survey.  That  survey 
was  taken  during  the  transition  from  a  relaxation  period  to  an  upwelling  pe¬ 
riod  and  is  therefore  much  more  similar  to  the  second  relaxation  state  than  to 
the  intervening  upwelling  state.  Another  factor  is  the  impact  of  open  boundary 
uncertainties  which  become  significant  almost  3  weeks  after  the  first  larger- 
scale  initialization  survey  and  just  before  the  second  re-initialization. 

Figures  10  and  11  show  the  temperature  RMS  differences  at  all  the  computed 
levels  for  the  one  and  two  day  forecasts  respectively.  The  temperature  RMS 
shows  greater  skill  and  variability  near  the  surface  than  at  depth.  A  similar 
trend  is  found  in  salinity  (not  shown).  At  30m,  the  skill  is  increased  after  Aug 
15  up  to  Aug  20  due  to  the  updates  in  model  parameters  (see  Lermusiaux, 
2007).  After  that,  the  skill  remains  close  to  zero  or  is  negative.  This  appears 
to  be  due  to  too  coarse  (3  to  5m)  vertical  grid  spacing  offshore,  near  the  bot¬ 
tom  of  the  summer  thermocline,  to  limited  data  assimilation  at  these  depths 
because  of  issues  with  some  of  the  real-time  salinity  data  and  to  tidal  effects 
which  are  neither  part  of  the  forecast  nor  the  persistence.  At  deeper  depths, 
the  decreased  variability  and  very  limited  data  leads  to  skill  estimates  less 
accurate  statistically.  An  interesting  result  is  that  forecast  skill  at  150  m  and 
300  m  remains  positive  about  one  week  before  the  second  initialization  survey 
is  available,  while  it  does  not  at  30  m.  This  may  be  because  below  the  thermo¬ 
cline,  boundary  uncertainties  have  not  yet  affected  the  larger-scale  fields.  The 
volume  mean  RMS  differences  most  closely  resembles  the  near  surface  differ- 
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ences.  This  reflects  both  the  larger  variations  at  the  surface  and  the  decreased 
number  of  data  values  at  depth. 


The  second  skill  metric  is  the  Pattern  Correlation  Coefficient  (PCC)  or  Anomaly 
Correlation  Coefficient.  Using  the  notation  for  the  RMS  errors,  the  forecast 
PCC  for  temperature,  TpeeForecast,  is  given  by 


TpccF, 


brecast 


(Tf  -  T„)  (T„  -  T„ 


Tf  —  Tr 


(1) 


In  Eq.  1  the  large-scale  mean,  T0,  is  removed  from  every  field.  Hence,  the  PCC 
measures  forecast  skill  based  on  the  correlations  of  anomalies  from  this  large- 
scale  mean.  In  this  work,  T0  is  a  weighted  mean  of  the  mapped  observations, 
with  the  weighting  function  e“(r/45)  where  r  is  the  distance  (km)  between  an 
observation  point  and  the  point  where  the  mean  is  being  evaluated.  This  mean 
is  computed  each  day,  using  only  that  day’s  mapped  observations.  The  scaling 
in  Eq.  1  is  such  that  TpeeForecast  £  [—1  1]  with  perfect  correlation  given  by 
TpeeForecast  =  1,  a  complete  lack  of  correlation  given  by  TpeeForecast  =  0  and 
perfect  anti-correlation  given  by  T  pec  Forecast  =  — 1-  As  in  the  RMS  error,  the 
persistence  PCC,  TpccpersiStence,  is  obtained  by  replacing  Tf  with  T;  in  Eq.  1. 
Using  PCC,  a  forecast  is  said  to  have  skill  if  Tpeeporecast  >  Tpccpersistence- 

Figure  12  shows  the  mean  PCC  differences  between  the  forecast  and  persis¬ 
tence  for  the  issued  real-time  one  and  two  day  forecasts.  As  with  the  RMS 
errors,  both  temperature  and  salinity  forecasts  have,  on  average,  skill  out  to 
2  days.  Here,  the  mean  temperature  PCC  differences  show  similar  skill  at  one 
and  two  days  (~0.16  mean  difference),  while  salinity  shows  greater  skill  at  2 
days  than  1  (0.072  mean  difference  at  one  day,  0.14  at  two  days).  The  depth 
dependence  (not  shown)  is  similar  to  the  RMS  error  dependence,  except  for 
the  fact  that  since  PCC  is  already  normalized,  the  PCC  amplitudes  do  not 
decrease  with  depth.  On  average,  PCC  skills  are  better  than  RMS  skills,  which 
indicates  that  the  simulation  captures  the  mesoscale  features  relatively  well. 


4  Reanalysis  Methodology 


The  real-time  forecasts  present  a  series  of  realizations  of  the  ocean  in  and 
around  Monterey  Bay.  These  realizations  are,  however,  somewhat  disjoint. 
This  arises  from  three  main  sources:  (1)  different  configurations  (stand-alone 
01,  nested-OI  and  stand-alone  ESSE)  were  employed  for  different  forecasts 
(figure  6),  (2)  different  starting  and  boundary  conditions  were  employed  at 
different  times  in  the  experiment  (historical  data  and  data  based  on  the  first 
and  second  R/V  Pt.  Sur  cruises),  and  (3)  model  parameters  were  being  tuned 
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during  the  experiment.  These  differences  prevent  the  real-time  forecasts  from 
being  simply  concatenated  together  and  used  for  serious  dynamical  studies. 
In  addition,  quality  control  of  the  various  data  sets  was  completed  the  year 
after  the  real-time  exercise.  Therefore  a  number  of  reanalyses  of  the  AOSN-II 
data  and  forecasts  were  performed  to  provide  a  single,  continuous  realization 
of  the  experiment.  We  next  present  the  methodology  of  the  final  re- analysis. 


The  goal  of  the  re-analysis  is  to  produce  a  better  tuned  continuous  4D  repre¬ 
sentation  of  the  oceans  in  the  AOSN-II  region  with  quality  controlled  data  as 
well  as  numerical  and  dynamical  parameter  improvements.  Because  of  the  lack 
of  offshore  larger-scale  data,  the  simplest  configuration  (stand-alone  Data  do¬ 
main  with  01  assimilation)  is  used  first.  To  create  the  initial  and  open  bound¬ 
ary  conditions,  three  objective  analyses  were  made  using  the  larger  mesoscale 
“initialization”  parameters  in  table  1.  The  first  was  an  analysis  for  OOZ  on  6 
Aug  2003  using  R/V  Pt.  Sur  CTDs,  WHOI  &  SIO  glider  and  NPS  SST  data 
in  the  period  2-6  Aug  2003.  The  second  was  an  analysis  for  OOZ  on  23  Aug 
2003  using  R/V  Pt.  Sur  and  R/V  John  Martin  CTDs,  WHOI  &  SIO  glider 
and  NPS  SST  data  in  the  period  21-25  Aug  2003.  The  third  was  an  analysis 
for  OOZ  on  5  Sep  2003  using  R/V  Pt.  Sur  CTDs  in  the  period  3-6  Sep  2003. 
These  analyses  were  then  processed  (with  a  level  of  no  motion  set  to  1250m) 
to  produce  initial  condition  fields  and  time-evolving  boundary  fields. 


The  boundary  fields  were  used  in  the  following  manner.  For  temperature  and 
salinity,  the  simulation  boundary  conditions  were  simply  linear  temporal  in¬ 
terpolations  of  the  analyses  at  the  boundary,  with  persistence  after  5  Sep. 
For  internal  mode  velocity,  the  time  interpolated  boundary  data  were  used 
as  the  “larger-scale  fields”  in  the  boundary  condition  relaxation  scheme  (ap¬ 
pendix  II).  The  transport  streamfunction  and  barotropic  vorticity  used  Or- 
lanski  radiation  conditions. 


Thirty-one  assimilation  fields  were  created  from  daily  objective  analyses  of  the 
August  6  -  September  06  2003  hydrographic  data  collected  by  the  R/Vs  Pt. 
Sur,  John  Martin  and  Pt.  Lobos;  the  WHOI  and  SIO  gliders;  and,  the  NPS 
aircraft  SST.  The  correlation  scales  were  those  of  table  1  and  analysis  dates  of 
7  Aug  -  6  Sep,  2003  at  1200Z.  For  each  analysis,  only  data  within  ±18  hours 
of  the  analysis  time  were  used.  From  these  mapped  temperature  and  salinity 
fields,  velocity  first-guess  estimates  were  obtained  by  geostrophic  integration. 
Only  the  temperature,  salinity  and  internal  velocities  were  assimilated,  follow¬ 
ing  the  schedule  and  weights  of  table  2,  with  the  provision  that  the  weights  for 
assimilating  internal  velocity  were  further  reduced  to  half  those  in  the  table. 
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5  Reanalysis  Results 


The  main  numerical  result  of  the  re-analysis  is  the  improvement  of  the  long¬ 
term  robustness  of  the  simulation,  primarily  by  controlling  instabilities  in 
the  southern  boundary  near  Point  Sur.  The  three  most  important  changes 
in  this  regard  were:  (1)  the  use  of  time-varying  boundary  data,  described  in 
the  section  4;  (2)  additional  smoothing  of  the  topography  in  the  neighborhood 
of  Point  Sur,  restricting  the  maximum  slope  to  0.11;  and  (3)  weakening  the 
relaxation  parameters  (table  3)  used  in  the  boundary  condition  relaxation 
scheme.  Additionally,  to  improve  coastal  currents,  the  coastal  friction  was 
weakened  by  increasing  the  relaxation  time-scale  (7200  s).  Otherwise,  the  run 
parameters  were  the  same  as  those  used  at  the  end  of  the  real-time  experiment. 

Figures  13  and  14  display  the  reanalysis  temperature  and  velocity,  at  the 
surface  and  30  nr  respectively,  every  2  days  in  the  period  10  Aug  -  1  Sep 
2003.  The  period  10-16  Aug  is  in  a  strong  upwelling  state.  Both  levels  show 
a  southward  flowing  coastal  current  which  crosses  the  mouth  of  Monterey 
Bay  and  is  deflected  west  (offshore)  north  of  Point  Sur.  The  surface  fields 
show  upwelling  plumes  from  Points  Ano  Nuevo  and  Sur.  The  plume  from  Ano 
Nuevo  extends  south  across  the  mouth  of  Monterey  Bay  and  merges  with  the 
plume  from  Point  Sur  by  Aug  13.  At  30  m,  a  broad  flow  enters  the  domain 
through  the  western  boundary  around  36°  12’N  and  turns  northwest.  Along  the 
northern  shelf/slope,  an  elongated  anticyclone  spans  the  region  between  the 
southward  coastal  current  and  the  northward  offshore  flow.  This  anticyclone 
retreats  to  the  north  while  the  offshore  branch  of  the  coastal  current  curves 
back  and  spins  up  another  anticyclone,  this  time  in  the  mouth  of  Monterey 
Bay.  Overall,  the  circulation  in  the  upper  layers  of  the  Bay  is  cyclonic  during 
such  upwelling  conditions. 

The  period  18-20  Aug  is  a  transition  between  upwelling  and  relaxation  states. 
The  winds  at  Ml  are  generally  weaker  and  their  direction  is  more  variable 
than  before  (figure  2).  The  upwelling  at  Ano  Nuevo  shuts  down  and  the  buoy¬ 
ancy  flow  patterns  of  the  thermocline  (30m)  begin  to  appear  also  at  the  sur¬ 
face.  At  30  m  the  two  previously  described  eddies  propagate  northward  (the 
northern  eddy  passing  through  the  domain  entirely).  The  period  22-26  Aug 
is  largely  a  relaxation  state  (with  the  26th  arguably  more  a  transition  day). 
The  surface  flows  remain  controlled  by  the  thermocline  flows.  During  this  pe¬ 
riod  in  the  real-time  experiment,  an  apparent  release  of  kinetic  energy  from 
the  wind-driven  surface  mean  circulation  to  internal  mesoscale  ocean  features 
(including  eddies  and  jets)  was  visualized  and  described  by  the  ESSE  uncer¬ 
tainty  error  fields  (Lermusiaux,  2006).  This  mechanism  was  confirmed  by  the 
detailed  multi-scale  energy  and  vorticity  dynamical  analysis  of  the  reanalysis 
fields  by  Liang  and  Robinson  (2008). 
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During  28-30  Aug,  the  winds  are  stronger  again  and  in  an  upwelling  con¬ 
figuration.  The  winds  drive  the  surface  flows  and  restart  upwelling  at  Point 
Ano  Nuevo  and  Point  Sur.  By  1  Sep,  winds  die  back  down  and  the  thermo- 
cline  flows  reassert  control  of  the  surface  flows.  An  interesting  characteristic 
to  notice  is  that  the  upwelling  patterns  of  the  first  and  second  upwelling  are 
different.  The  first  upwelling  lead  to  wide  westward  plumes  at  the  two  Points, 
while  the  second  does  not.  During  28-30  Aug,  the  upwelling  front  is  relatively 
parallel  to  the  coast,  just  west  of  the  mouth  of  Monterey  Bay. 

To  assess  the  skill  of  the  re-analysis,  the  simulation  was  repeated  with  the 
assimilation  of  data  stopping  on  20  August.  Skill  metrics,  as  described  in  the 
real-time  section,  were  then  constructed  comparing  the  repeat  simulation  to 
the  fields  that  would  have  been  assimilated  on  21-25  August  (corresponding 
to  the  second  R/V  Pt.  Sur  survey).  For  these  comparisons,  the  persistence 
fields  were  constructed  from  a  time  dependent  objective  analysis  (OA)  of  all 
data  prior  to  21  August  with  an  analysis  date  of  20  August  (he.  a  “best, 
data  only”  estimate  for  20  August).  This  persistence  OA  was  made  using  the 
(larger)  initialization  scales  (table  1). 

Figure  15  shows  the  resulting  differences  in  RMS  errors  between  the  temper¬ 
ature  fields  of  the  persistence  and  the  above-described  re-analysis  simulation 
without  data  assimilation  after  20  August.  For  this  reanalysis,  the  skill  is 
positive  at  all  depths  and  times,  except  on  3  instances.  This  is  a  large  im¬ 
provement  when  compared  to  the  real-time  forecast.  It  is  the  quality  control 
of  the  ocean  data  and  the  improved  model  parameters,  initial  and  boundary 
conditions  that  lead  to  this  large  improvement.  The  largest  skill  in  amplitude 
mainly  resides  in  the  thcrmocline  (30m),  which  is  the  location  of  the  largest 
dynamical  variability  during  relaxation  conditions.  The  skill  at  the  surface  is 
smaller  in  amplitude  in  part  because  of  this.  It  is  also  smaller  than  for  the 
real-time  forecast  because  of  the  closer  temporal  proximity  of  the  20  Aug  field 
persistence  fields  to  the  new  21-25  Aug  data. 

A  refined  look  at  the  vertical  error  structure  was  made  by  comparing  the  R/V 
Pt.  Sur  CTDs  from  the  second  survey  to  model-simulation  profiles  obtained 
by  nearest-neighbor  interpolation  (space  and  time  from  Aug  21  to  25)  from 
the  persistence  fields  and  the  repeat  re-analysis  simulation  fields  (with  no  DA 
after  Aug  20).  These  profiles  were  then  vertically  interpolated  to  a  common 
1  m  vertical  grid  and  differenced  (simulation  -  data;  persistence  -  data).  These 
error  profiles  are  then  averaged  over  all  profiles  using  either  a  simple  mean 
(bias)  or  an  root- mean-square  estimate.  The  resulting  mean  temperature  error 
profiles  are  presented  in  figure  16.  Again,  skill  is  seen  in  the  main  thermocline 
down  to  roughly  the  bottom  of  the  California  Under  Current  (300m).  At  all  of 
these  dynamically  active  depths,  the  persistence  “data-only”  estimate  is  not 
as  good  as  the  model  prediction.  Below  300m,  there  is  a  nearly  linear  increase 
in  the  bias  down  to  500m,  followed  by  a  more  gradual  decrease.  This  loss  of 
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skill  is  clue  in  part  to  the  general  lack  of  data  below  400m  between  the  R/V 
Pt.  Snr  surveys  (no  glider  data  below  400m  and  only  intermittent  small  CTD 
surveys).  Even  though  there  is  very  little  data  to  understand  the  cause  of  this 
deep  drift  of  the  model  estimates,  it  could  possibly  be  corrected  by  a  better 
tuning  of  the  deep  model  parameters. 


6  Conclusion 


The  methodologies  and  results  of  the  utilization  of  HOPS  during  AOSN-11 
were  presented  and  studied.  These  included  a  description  of  the  modeling 
approach  and  modeling  system  components.  The  approach  consisted  of  ini¬ 
tializing  the  modeling  system  from  adequate  background  ocean  fields,  which 
is  often  obtained  from  a  synoptic  survey  of  the  larger-scale  conditions  in  the 
region.  Subsequently,  data  are  assimilated,  the  model  skill  is  evaluated  and 
model  properties  are  improved,  all  in  real-time.  Importantly,  the  model  out¬ 
puts  are  studied  in  as  much  details  as  possible  within  the  time  and  personnel 
constraints,  so  as  to  foster  real-time  improvements  in  possibly  all  components 
of  the  modeling  system. 

During  AOSN-ff,  conditions  were  on  average  in  anomalously  strong  upwclling 
favorable  conditions.  Two  week-to-10-days  long  upwelling  events  were  sepa¬ 
rated  by  two  3-to-5-days  long  relaxation  events.  HOPS  was  re-initialized  three 
times,  once  from  historical  data,  twice  from  each  of  the  first  two  R/V  Pt. 
Sur  surveys.  The  modeling  system  assimilated  various  data  types  (Ship  CTD, 
AUVs,  gliders  and  SST  data)  and  model  parameters  were  improved  in  real¬ 
time.  From  the  real-time  results,  one  can  already  obtain  a  description  of  several 
general  circulation  and  hydrographic  properties  of  the  region  that  remained 
valid  after  re-analysis  and  model  improvements. 

The  qualitative  evaluation  of  the  real-time  forecast  skill  showed  that  many 
of  the  surface  ocean  features  were  predicted,  but  that  their  detailed  position 
and  shapes  as  well  as  their  (sub)-mesoscale  properties  were  often  not  accu¬ 
rate,  especially  offshore  where  data  were  limited.  In  the  Bay,  the  match  with 
surface  velocity  data  (not  assimilated)  was  on  average  good.  This  indicates 
that  the  WHOI  gliders  measuring  the  northern  portion  of  the  Bay  allowed  the 
estimation  of  the  surface  circulation  over  the  whole  Bay,  even  in  the  absence  of 
tidal  forcing.  This  was  feasible  because  of  the  use  of  a  data-assimilative  ocean 
model,  including  the  other  data  sets  (see  Fig.  3)  and  the  atmospheric  predic¬ 
tions.  Note  that  the  HOPS  surface  current  predictions  in  the  Bay  were  also 
found  accurate  enough  to  estimate  the  LCS  in  a  very  useful  way  for  surface 
drifter  planning  and  tracking. 

The  quantitative  evaluation  of  the  real-time  forecast  skill  based  on  RMS  error 
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comparisons  showed  that  the  real-time  forecasts  had,  on  average,  skill  out  to 
2  days.  The  skill  was  found  to  be  greater  near  the  surface.  The  PCC  com¬ 
parisons  showed,  on  average,  greater  skill  than  the  RMS  comparisons.  This 
indicates  that  the  mesoscale  features  were  overall  well  represented  in  the  real¬ 
time  forecast.  For  simulations  of  10  days  or  more  in  duration,  it  was  found 
that  uncertainties  in  the  open  boundary  forcing  of  our  Data  Domain  could 
lead  to  significant  errors  in  the  Monterey  Bay  region. 

The  main  results  of  the  reanalysis  effort  is  an  improvement  of  the  long-term 
stability  of  the  simulation  through  improved  model  parameters  and  open 
boundary  conditions.  These  boundary  conditions  are  based  on  a  time- varying 
larger-scale  objective  analysis  of  observations.  The  quantitative  skill  of  the 
re-analysis  was  also  measured  and  improved  when  compared  to  the  real-time 
fields.  These  reanalysis  fields  allowed  a  description  of  the  significant  ocean 
features  that  occurred  during  AOSN-II.  Two  types  of  upwelling  events  were 
observed,  one  with  bifurcations  or  plumes  extending  westward  at  Point  AN 
and  Point  Sur,  and  one  with  a  thinner  band  of  upwelled  water  parallel  to  the 
coast  and  across  Monterey  Bay.  During  strong  upwelling  events  the  flow  in 
the  upper  10-20  meters  has  scales  similar  to  atmospheric  scales,  with  broad 
features.  Once  the  winds  subside,  kinetic  energy  becomes  available  and  leads 
to  mesoscale  features  developing  within  a  warming  (relaxing)  upper  thermo- 
cline.  The  estimated  flow  at  deeper  100  to  300m  depths  consisted  on  average 
of  an  overall  broad  northward  flow,  sometimes  with  a  branch  following  the 
coastal  topographic  features.  Also  simulated  were  the  northward  motions  of 
the  anticyclone  and  other  eddies  that  often  develop  in  front  of  Monterey  Bay, 
and  the  westward  plumes  that  develop  in  upwelling  conditions. 

A  number  of  important  lessons  for  future  work  were  learned  from  this  exper¬ 
iment.  The  first  lesson  was  based  on  the  existence  of  a  strong  daily  cycle  in 
the  upper  layers  and  on  the  shelf.  To  improve  modeling  this  cycle,  the  data 
need  to  be  assimilated  more  frequently,  with  shorter  time  groupings  (e.g.  twice 
daily  assimilations,  at  the  times  of  maximum/minimum  observed  SST,  with 
correspondingly  smaller  time  windows,  §4  &  table  2,  and  correlation  scales, 
table  1).  A  second  lesson  was  the  importance  of  sufficient  synoptic  accuracy 
in  the  larger  domain  for  two-way  nesting.  Feature  models  could  be  used  to  ex¬ 
tend  limited  observations  in  a  dynamically  useful  way  for  initialization  of  the 
larger  domain.  Output  from  another  dynamical  model  (with  sufficient  accu¬ 
racy)  could  be  another  source  for  initial/boundary  conditions.  A  third  lesson 
dealt  with  the  specifics  of  the  “Data  Domain” .  A  finer  resolution  near  the  ther- 
mocline  was  indicated  by  the  evaluation  of  the  real-time  forecasts  (§3.3).  The 
prevalence  of  numerical  issues  at  the  southern  boundary  argued  for  moving 
south,  to  include  the  entire  dynamic  upwelling  center  region  in  the  numeri¬ 
cal  domain.  A  fourth  lesson  came  from  the  results  of  the  reanalysis,  namely  a 
continuous  set  of  fields  covering  the  experiment,  and  the  improved  skill  that  re¬ 
sulted.  This  could  be  realized  in  a  real-time  experiment  by  performing  longer, 
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full  experiment  sensitivity  studies  overnight  (between  the  issuance  of  the  prod¬ 
ucts  and  the  start  of  the  next  day’s  operations,  Fig.  4).  There  is  additional 
work  and  research  remaining  in  studies  of  the  AOSN-II  data  based  on  dynam¬ 
ical  models.  One  should  also  investigate  the  effects  of  adding  a  hierarchy  of 
ever  more  complex  tidal  parameterizations,  from  a  simple  increase  of  mixing 
coefficients  as  a  function  of  barotropic  tidal  amplitudes  to  a  full  free  surface 
simulations  forced  with  barotropic  tides  (Rosenfeld  et  ah,  2008;  Wang  et  al., 
2008).  Some  of  the  data  sets  which  were  not  assimilated  in  the  re-analysis 
presented  could  be  used;  they  include  buoy,  CODAR  (Shulman  and  Paduan, 
2008)  and  SSH  data.  Learning  from  AOSN-II,  most  of  these  updates  were 
carried  out  in  August  2006  for  the  real-time  Monterey  Bay  06  experiment. 
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7  Appendix  I:  Characteristics  of  the  Real-time  Modeling  System 
Components  and  Processing 


Detailed  characteristics  of  the  components  and  processing  of  the  real-time 
modeling  system  are  described  next.  These  include  the  modeling  domains, 
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data  processing  and  analyses,  atmospheric  forcing  and  product  dissemina¬ 
tion.  Parameterizations,  parameters  and  data  assimilation  are  discussed  in 
the  Sect.  2.3. 


7.1  Domains,  Nesting,  Grids  and  Bathymetry 


The  Data  Domain  was  centered  at  (36.62N,  122. 40W)  with  an  83x96  grid  at 
1.5km  resolution.  The  Off  Shore  Domain  was  centered  at  (36.18N,  123. 59W) 
with  an  83x96  grid  at  4.5km  resolution.  Both  domains  had  22  terrain- following 
vertical  levels  arranged  in  a  double-sigma  configuration  (Haley,  1999)  with  12 
levels  in  the  upper  system  and  10  in  the  lower.  Two-way  interprocess  commu¬ 
nications  exchange  information  between  the  two  model  domains.  Specifically, 
the  Off  Shore  Domain  in  the  region  of  overlap  is  replaced  with  information 
averaged  from  the  Data  Domain  and  the  boundaries  of  the  Data  Domain  are 
constrained  by  information  interpolated  from  the  updated  Off  Shore  Domain. 

To  obtain  a  bathymetry  for  the  domains,  a  30  arc  second  topography  was  con¬ 
structed  a  follows.  First  a  25m  gridded  topography9  (MBARI,  2000),  covering 
36.44-37.06N  122. 51-121. 78W  was  averaged  to  a  30  arc  second  grid.  Outside 
this  grid,  elevations  were  obtained  from  the  30  arc  second  GLOBE  10  dataset 
(Hastings  and  Dunbar,  1999).  Bathymetry  outside  the  range  of  the  MBARI 
dataset,  but  inside  35.52-37.99N  123. 975-121. 041W  were  obtained  from  a  0.01° 
dataset  constructed  by  NPS  11 .  Bathymetry  outside  of  the  range  of  the  NPS 
dataset  were  taken  from  a  1  minute  DBDB-V 12  dataset  (Steed  et  ah,  2002). 

The  resultant  topography  was  then  bilinearly  interpolated  to  the  horizontal 
grid  of  each  domain.  The  interpolated  topographies  were  then  conditioned 
by:  (1)  resetting  all  depths  shallower  than  10m  to  10m;  (2)  applying  a  two 
dimensional  median  filter  (Lozano  et  al.,  1994);  and  (3)  applying  a  smoothing 
algorithm  to  simultaneously  bound  the  reduced  slope  by  2.5  and  the  slope  by 
0.225  (Haley  and  Lozano,  2001). 


9  provided  by  John  Ryan 

1(1  http :  //www.ngdc .  noaa.gov/mgg/topo/globe .  html 

11  “Derived  from  a  combination  of  two  sources:  (1)  Monterey  Bay  Aquarium  Re¬ 
search  Institute,  Moss  Landing,  California  and  (2)  U.S.  Geological  Survey,  Menlo 
Park,  California”  Jeff  Paduan,  NPS  (personal  communication) 

12  https : // 128 . 160 . 23 . 42/dbdbv/dbvquery .html 
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7.2  Oceanographic  Data  Processing  and  Objective  Analyses 


On  a  daily  basis,  the  available  CTD  profiles  and  WHOI  and  SIO  glider  pseudo¬ 
profiles  were  collected  and  converted  to  the  HOPS  ASCII  format.  First,  the 
data  from  each  individual  instrument  were  manually  examined  for  internal 
consistency.  When  needed,  a  median  filter  (Tukey,  1977)  was  applied  to  de¬ 
spike  noisy  data.  Then  the  profiles  from  different  instruments  were  manually 
compared  to  each  other  and  to  the  most  recent  R/V  Pt.  Sur  CTD  survey. 
Inconsistent  data  points  were  identified  and  removed.  Every  effort  was  made 
to  retain  as  much  data  as  possible  (e.g.  If  a  profile  fell  within  the  common 
envelope  in  its  surface  and  deep  values  but  was  far  off  in  the  thermocline,  only 
the  thermocline  portion  would  be  removed.). 

The  in  situ  and  SST  profiles  were  gathered  together  and  processed  through 
the  HOPS  objective  analysis  package.  The  data  were  then  analyzed  with  the 
correlation  parameters  in  table  1,  depending  on  whether  the  analyzed  fields 
were  to  be  used  to  initialize  a  simulation  or  to  be  assimilated  into  a  simulation. 
The  outputs  of  the  objective  analysis  were  three  dimensional  fields  of  temper¬ 
ature,  salinity,  dynamic  height  and  their  associated  errors.  These  fields  were 
located  on  the  same  horizontal  grid  as  the  PE  domain,  but  were  on  flat  lev¬ 
els  in  the  vertical.  The  temperature  and  salinity  were  interpolated  to  the  PE 
terrain-following  vertical  coordinates.  The  dynamic  height  was  first  differenti¬ 
ated  to  produce  horizontal  velocities  on  the  flat  levels.  These  velocities  were 
then  interpolated  to  the  terrain-following  coordinates,  and  decomposed  into 
a  depth  averaged  component  and  an  “internal”  (remainder)  component.  The 
depth  averaged  component  was  used  to  generate  a  transport  streamfunction, 
to  ensure  the  non- divergence  of  transport. 


7.3  Atmospheric  Forcing  Flux  Processing 

The  Navy’s  operational  Coupled  Ocean/Atmosphere  Mesoscale  Prediction 
System  (COAMPS)  analysis  and  72  hour  forecast  fields  (Doyle  et  ah,  2008) 
were  processed  prior  to  their  utilization.  The  wind  stress  fields  from  COAMPS 
were  bi-cubically  interpolated  to  the  HOPS  grids,  rescaled  from  (N/m2)  to 
(dyne/cm2)  and  decomposed  in  terms  of  the  direction  vectors  aligned  with 
the  HOPS  grids.  Special  care  was  taken  to  resolve  the  mismatch  of  HOPS  and 
COAMPS  land  masks.  COAMPS  wind  stress  values  are,  on  average,  an  order 
of  magnitude  higher  over  the  land  than  over  the  sea.  Hence,  discrepancies  in 
the  land  masks  could  lead  to  excessive  values  of  the  wind  stress  on  the  HOPS 
grids  near  the  coasts.  To  control  this  problem,  only  wind  stress  values  from 
COAMPS  sea-points  were  used  to  interpolate/extrapolate  to  the  HOPS  grids. 
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The  surface  net  heat  flux  into  the  ocean  (W/m2)  was  composed  as  the  sum 
of  the  shortwave  radiation,  the  longwave  radiation,  the  sensible  heat  flux  and 
the  latent  heat  flux;  all  of  which  were  directly  available  from  the  COAMPS 
output  and  bi-cubically  interpolated  to  the  HOPS  grids. 

The  net  surface  water  flux  out  of  the  ocean  (cm/day)  was  calculated  as 
the  difference  of  evaporation  minus  precipitation.  Evaporation  was  obtained 
by  dividing  the  COAMPS  latent  heat  flux  by  the  latent  heat  of  evapora¬ 
tion.  The  latent  heat  of  evaporation,  Lv  (J/kg),  was  computed  from  the 
COAMPS  air  temperature,  Tair  (°  C),  using  the  empirical  formula  Lv  = 
2.5008(106)  —  2.3(103)Tair  (Gill,  1982,  p607).  The  resulting  evaporation  was 
then  converted  from  (kg/(m2s))  to  (cm/day).  The  total  precipitation  was  ob¬ 
tained  directly  from  COAMPS,  merely  needing  a  conversion  from  (mm  in  12 
hours)  to  (cm/day).  The  net  water  flux  was  bi-cubically  interpolated  to  the 
HOPS  grids. 

The  validation  of  COAMPS  re-analysis  products  for  the  Monterey  Bay  was 
examined  by  Kindle  et  al.  (2002)  using  the  Monterey  Bay  Aquarium  Research 
Institute  (MBARI)  mooring  wind  measurements.  A  two-year  time  series  of 
observed  winds  at  MBARI  moorings  Ml  and  M2  were  found  to  be  in  good 
agreement  with  the  model  outputs,  both  on  the  synoptic  scale  and  for  summer¬ 
time  diurnal  variability.  During  AOSN-II,  we  compared  a  two-week  sequence  of 
COAMPS  wind  forecasts  (with  lead  times  of  24-hours,  48-hours  and  72-hours) 
to  the  observational  data  from  MBARI  moorings  Ml  and  M2.  The  forecasts 
at  all  examined  lead  times  had  winds  with  a  generally  correct  magnitude 
(±2m/s),  reproduced  well  the  observed  diurnal  cycle  at  both  Ml  and  M2,  and 
consistently  reproduced  the  observed  strong  wind  and  wind  relaxation  events. 


7.4  Product  Dissemination 


Plots  of  the  HOPS  real-time  nowcast /forecast  fields  were  uploaded  to  a  local 
web  page  (§2.2).  The  horizontal  plots  were  provided  at  4  depths,  pre-selected 
based  on  the  vertical  structure  of  the  historical  data:  surface,  10m  (base  of 
mixed  layer),  30m  (middle  of  thermocline),  and  200m  (California  Under  Cur¬ 
rent  region).  Vertical  sections  were  made  coming  out  of  points  A110  Nuevo  and 
Sur  and  out  of  the  center  of  Monterey  Bay.  Additional  links  include  special 
products,  pre-experiment  planning  documents  and  post-experiment  on-going 
analyses. 

NetCDF  output  files  from  the  issued  HOPS  PE  forecasts  were  uploaded  to 
the  MBARI  server  (§2.2).  Specially  requested  products  were  also  served  there 
(e.g.  output  for  LCS  analysis). 
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8  Appendix  II:  New  Open  Boundary  Conditions 


The  southern  open  boundary  by  Point  Sur  was  especially  problematic  and 
limiting  to  simulation  duration.  To  handle  this  boundary,  a  variation  on  the 
standard  Orlanski  radiation  condition  was  implemented  and  used  in  real-time 
(see  also  Lermusiaux,  2007).  It  is  based  on  relaxing  the  predicted  boundary 
values  to  a  control  field  and  is  described  next.  Note  that  a  set  of  open-boundary 
conditions  that  are  variations  of  this  Orlanski-relaxation  theme  have  also  been 
implemented  and  will  be  reported  elsewhere. 


8.1  Orlanski  and  relaxatio7i  based  open  boundary  conditions 


“Boundary  condition  relaxation”  is  a  relaxation-type  forcing  applied  to  the 
predicted  open  boundary  values  (e.g.  by  an  Orlanski  scheme).  This  forcing 
imposes  a  tunable  and  weak  constraint  towards  some  background  conditions. 
Of  course,  these  background  conditions  can  change  with  time,  for  example  due 
to  a  larger-scale  buoyancy  forcing.  The  basic  form  of  this  forcing  is 

=  ubcbXl  ~  ^b,k  -  udatabi  k)  (2) 


where  u  is  the  variable  whose  boundary  values  are  being  sought,  b  is  the 
boundary  index,  i  is  the  along-boundary  index,  k  is  the  depth  index,  ubc 
is  the  predicted  boundary  value,  Udata  is  some  additional  data  for  u  along 
the  boundaries  and  l ob^  is  the  relaxation  weight.  The  relaxation  weight  is 
computed  from 
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Tlb  +  t2  b 
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where  At  is  the  time  step  (s),  T\  is  the  base  decay  rate  (s),  r2  is  the  depth  decay 
rate  (s),  7  is  a  nondimensional  depth  decay  factor  and  K  is  the  total  number 
of  model  levels.  For  model  levels  shallower  than  K  —  kfiX,  the  local  velocity 
ubil  is  examined  and,  if  found  to  be  outflow,  the  weight  u>b,k  is  multiplied  by 
All  of  these  parameters  were  tuned  to  onr  specific  situation  in  AOSN-II  and 
we  expect  that  their  values  should  be  modified  for  other  periods  or  regions. 


References 

Bellingham,  J.G.  and  Y.  Zhang  (2005)  Observing  processes  that  vary  in  time 
and  space  with  heterogeneous  mobile  networks.  Proceedings  of  the  Interna- 


25 


tional  Workshop  on  Underwater  Robotics,  Genoa,  Italy,  pp.  9-16. 

Bellingham,  J.G.  (2006)  New  challenges  for  ocean  technology.  Marine  Tech¬ 
nology  Society  Journal,  40.  139-140. 

Bretherton,  F.P.,  R.E.  Davis  and  C.B.  Fandry  (1976)  A  Technique  for  Objec¬ 
tive  Analysis  and  Design  of  Oceanographic  Experiments  Applied  to  MODE- 
73.  Deep  Sea  Research,  23.  559-582. 

Brink,  K.H.,  R.C.  Beardsley,  P.P.  Niiler,  M.  Abbot,  A.  Huyer,  S.  Ramp,  T. 
Stanton,  and  D.  Stuart.  (1991)  Statistical  properties  of  near-surface  flow  in 
the  California  coastal  transition  zone.  J.  Geophys.  Res.,  96.  14,693-14,706. 

Bryan,  K.  (1969)  A  Numerical  Method  for  the  Study  of  the  Circulation  of 
the  World  Ocean.  Journal  of  Computational  Physics,  4.  347-376. 

Carter,  E.F.  and  A.R.  Robinson  (1987)  Analysis  Models  for  the  Estimation 
of  Oceanic  Fields.  Journal  of  Atmospheric  and  Oceanic  Technology,  4(1). 
49-74. 

Chang,  K.  (2003)  Taking  the  Oceans’  Pulse,  With  Help  From  Robot  Subs. 
The  New  York  Times,  September  30,  2003.  Section  F;  Column  3;  Science 
Desk;  Pg.  1-2. 

Chao,  Y.,  Z.  Li,  J.  Farrara,  J.C.  McWilliams,  J.  Bellingham,  X.  Capet,  F. 
Chavez,  J.-K.  Choi,  R.  Davis,  J.  Doyle,  D.  Fratantoni,  P.  Li,  P.  Marchesiello, 
M.A.  Moline,  J.  Paduan  and  S.  Ramp  (2008)  Development,  Implementa¬ 
tion  and  Evaluation  of  a  Data-Assimilative  Ocean  Forecasting  System  off 
the  Central  California  Coast.  Deep  Sea  Research  II,  (submitted). 

Charney,  J.G.,  R.  Fjortoft,  and  J.  Von  Neumann  (1950)  Numerical  integra¬ 
tion  of  the  barotropic  vorticity  equation.  Tellus,  2.  237-254. 

Chavez,  F.P.,  R.T.  Barber,  P.M.  Kosro,  A.  Huyer,  S.R.  Ramp,  T.  Stanton  and 
B.  Rojas  de  Mendiola  (1991)  Horizontal  advection  and  the  distribution 
of  nutrients  in  the  coastal  transition  zone  off  northern  California:  Effects 
on  primary  production,  phytoplankton  biomass,  and  species  composition. 
Journal  of  Geophysical  Research ,  96,  14833-14848. 

Chavez,  F.P.,  J.T.  Pennington,  R.  Herlien,  H.  Jannasch,  G.  Thurmond,  and 
G.E.  Friederich  (1997)  Moorings  and  drifters  for  real-time  interdisciplinary 
oceanography,  Journal  of  Atmospheric  and  Ocean  Technology,  14.  1199- 
1211. 

Chavez,  F.P.  and  C.A.  Collins  (2000)  Studies  of  the  California  current  system: 
Present,  past  and  future.  Deep-Sea  Research  II,  47.  761-763. 

Collins,  C.  A.,  N.  Garfield,  R.G.  Paquette,  and  E.  Carter  (1996)  Lagrangian 
Measurement  of  subsurface  poleward  Flow  between  38°N  and  43°N  along 
the  West  Coast  of  the  LInited  States  during  Summer,  1993.  Geophys.  Res. 
Letters,  23(18).  2461-2464. 

Collins,  C.A.,  J.T.  Pennington,  C.G.  Castro,  T.A.  Rago  and  F.P.  Chavez 
(2003)  The  California  Current  system  off  Monterey,  California:  physical 
and  biological  coupling.  Deep  Sea  Res.  II,  50.  2389-2404. 

Curtin,  T.B.  and  J.G.  Bellingham  (2001)  Autonomous  ocean-sampling  net¬ 
works.  IEEE  Journal  of  Oceanic  Engineering,  26(4).  421-423. 

Curtin,  T.B.  and  J.G.  Bellingham  (2008)  Progress  Toward  Autonomous 


26 


Ocean  Sampling  Networks.  Deep  Sea  Research  II,  (submitted). 

Curtin,  T.B.,  J.G.  Bellingham,  J.  Catipovic,  and  D.  Webb  (1989)  Au¬ 
tonomous  oceanographic  sampling  networks.  Oceanography,  6.  86-94. 

Cushman- Roisin,  B.  (1994)  Introduction  to  geophysical  fluid  dynamics.  En¬ 
glewood  Cliffs,  NJ.  Prentice  Hall. 

Davis,  R.E.,  C.E.  Eriksen  and  C.P.  Jones  (2002)  Autonomous  buoyancy- 
driven  underwater  gliders.  In  The  Technology  and  Applications  of  Au¬ 
tonomous  Underwater  Vehicles.  (G.  Griffiths,  ed),  Taylor  and  Francis,  Lon¬ 
don.  324  pp. 

Davis,  R.E.,  N.  Leonard  and  D.M.  Fratantoni  (2008)  Routing  Strategies  for 
Underwater  Gliders.  Deep  Sea  Res.  II,  (submitted). 

Davis,  T.M.,  K.A.  Countryman  and  M.J.  Carron  (1986)  Tailored  acoustic 
products  utilizing  the  NAVOCEANO  GDEM  (a  generalized  digital  envi¬ 
ronmental  model),  in  Proceedings,  36th  Naval  Symposium  on  Uriderwater 
Acoustics,  Naval  Ocean  Systems  Center,  San  Diego,  CA. 

Doyle,  J.D.,  Q.  Jiang,  Y.  Chao,  and  J.  Farrara  (2008)  High-resolution  real¬ 
time  modeling  of  the  marine  atmospheric  boundary  layer  in  support  of  the 
AOSNII  field  campaign.  Deep  Sea  Research  II,  (submitted). 

Fiorelli,  E.,  N.E.  Leonard,  P.  Bhatta,  D.  Paley,  R.  Bachmayer,  and  D.M. 
Fratantoni  (2006)  Multi- AUV  control  and  adaptive  sampling  in  Monterey 
Bay.  IEEE  Journal  of  Oceanic  Engineering,  31(4).  935-948. 

Fox,  A.D.  and  S.J.  Masked  (1995)  Two-Way  Interactive  Nesting  of  Primitive 
Equation  Ocean  Models  with  Topography.  J.  Phys.  Oceanog.,  25.  2977- 
2996. 

Garfield,  N.,  C.A.  Collins,  R.G.  Paquette,  and  E.  Carter  (1999)  Lagrangian 
Exploration  of  the  California  Undercurrent,  1992-1995.  J.  Phys.  Oceanoqr., 

29.  560-583. 

Gill,  A.E.  (1982)  Atmosphere-Ocean  Dynamics.  Academic  Press,  New  York. 

662p. 

Haley,  P.J.  (1999)  Grids  User’s  Guide.  Harvard  Open  Ocean  Model  Reports, 

54,  Harvard  LIniversity,  Cambridge,  MA. 

Haley,  P.J.,  P.F.J.  Lermusiaux,  W.G.  Leslie  and  A.R.  Robin¬ 
son  (1999)  Harvard  Ocean  Prediction  System  (HOPS), 
http : / / oceans . deas . harvard . edu/HOPS/HOPS . html. 

Haley,  P.J.,  and  C.  J.  Lozano  (2001)  COND_TOPO:  Topography  conditioning 

in  Matlab.  http :  // oceans  .  deas .  harvard .  edu/HOPS/ Cond_Topo/ cond_topo .  ps .  gz 

Haller,  G.  (2002)  Lagrangian  coherent  structures  from  approximate  velocity 
data.  Phys.  Fluids,  14.  1851-1861. 

Hastings,  D.A.,  and  P.K.  Dunbar.  (1999)  Global  Land  One- km  Base  Elevation 
(GLOBE)  Digital  Elevation  Model  Documentation.  Publication  KGRD  3f, 

Boulder,  Colorado,  NOAA  National  Geophysical  Data  Center. 

Hayward,  T.L.  and  A.W.  Mantyla  (1990)  Physical,  chemical  and  biological 
structure  of  a  coastal  eddy  near  Cape  Mendocino.  J.  Mar.  Res.,  4.  825-850. 

Hickey,  B.M.  (1998)  Coastal  Oceanography  of  Western  North  America  from 
the  tip  of  Baja  California  to  Vancouver  Is.  In  The  Sea,  (K.H.  Brink  and 


27 


A.R.  Robinson,  eds.)  11.  Wiley  and  Sons,  Inc.  345-393. 

Hodnr,  R.  M.  (1997)  The  Naval  Research  Laboratory’s  Coupled 
Ocean/Atmosphere  Mesoscale  Prediction  System  (COAMPS).  Mon.  Wea. 
Rev.,  125,  1414-1430. 

Hnyer,  A.E.  (1983)  Coastal  upwelling  in  the  California  Current  system.  Progr. 
Oceanogr.  12.  259-284. 

Huyer,  A.,  J.A.  Barth,  P.M.  Kosro,  R.K.  Shearman,  and  R.L.  Smith  (1998) 
Llpper-ocean  water  mass  characteristics  of  the  California  Current,  Summer 
1993.  Deep  Sea  Res.  II,  45.  1411-1442. 

Huyer,  A.,  P.M.  Kosro,  J.  Flcischbein,  S.R.  Ramp,  T.  Stanton,  L.  Washburn, 
F.P.  Chavez,  T.J.  Cowles,  S.D.  Pierce,  and  R.L.  Smith  (1991)  Currents 
and  water  masses  of  the  Coastal  Transition  Zone  off  northern  California, 
June  to  August,  1988.  J.  Geophys.  Res.,  96.  14,809-14,831. 

Kindle,  J.C.,  R.M.  Hodur,  S.  deRada,  J.D.  Paduan,  L.K.  Rosenfcld,  and  F. 
Chavez  (2002)  A  COAMPS  reanalysis  for  the  Eastern  Pacific:  Properties 
of  the  diurnal  sea  breeze  along  the  central  California  coast.  Geophysical 
Research  Letters,  29(24).  56(1-4). 

Kosro,  P.M.,  A.  Huyer,  S.R.  Ramp,  R.L.  Smith,  F.P.  Chavez,  T.L.  Cowles, 
M.R.  Abbott,  P.T.  Strub,  R.T.  Barber,  P.  Jessen  and  L.F.  Small  (1991) 
The  structure  of  the  transition  zone  between  coastal  waters  and  the  open 
ocean  off  northern  California,  winter  and  spring  1987.  J.  Geophys.  Res.,  96. 
14,707-14,730. 

Lekien  F.,  Coulliette  C.,  Mariano  A.J.,  Ryan  E.H.,  Shay,  L.K.,  Haller,  G.  and 
Marsden  J.  (2005)  Pollution  Release  Tied  to  Invariant  Manifolds:  A  Case 
Study  for  the  Coast  of  Florida.  Physica  D,  210(1-2)  1-20. 

Leonard,  N.E.  and  J.G.  Graver  (2001)  Model-based  feedback  control  of  au¬ 
tonomous  underwater  gliders,  IEEE  Journal  of  Oceanic  Engineering,  26(4). 
633-645,  Special  Issue  on  Autonomous  Ocean-Sampling  Networks. 

Leonard,  N.E.,  D.  Paley,  F.  Lekien,  R.  Sepulchre,  and  D.M.  Fratantoni  (2007) 
Collective  motion,  sensor  networks  and  ocean  sampling.  Proceedings  of  the 
IEEE,  95(1)  48-77. 

Lermusiaux,  P.F.J.  (1997)  Error  Subspace  Data  Assimilation  Methods  for 
Ocean  Field  Estimation:  Theory,  Validation  and  Applications.  PhD  Thesis, 
Harvard  LIniversity,  Cambridge,  MA. 

Lermusiaux,  P.F.J.  (1999)  Data  assimilation  via  Error  Subspace  Statistical 
Estimation.  Part  II:  Middle  Atlantic  Bight  shelfbreak  front  simulations  and 
ESSE  validation.  Month.  Weather  Rev.,  127(7).  1408-1432. 

Lermusiaux  P.F.J.  (2001)  Evolving  the  subspace  of  the  three-dimensional 
multiscale  ocean  variability:  Massachusetts  Bay.  J.  Marine  Systems,  29(1- 
4).  385-422. 

Lermusiaux  P.F.J.  (2006)  Uncertainty  Estimation  and  Prediction  for  Inter¬ 
disciplinary  Ocean  Dynamics.  J.  Computational  Physics,  217(1).  176-199. 

Lermusiaux,  P.F.J.  (2007)  Adaptive  modeling,  adaptive  data  assimilation 
and  adaptive  sampling.  Physica  D,  230(1-2).  172-196. 

Lermusiaux  P.F.J.,  C.-S.  Chiu,  G.G.  Gawarkiewicz,  P.  Abbot,  A.R.  Robinson, 


R.N.  Miller,  P.J.  Haley,  W.G.  Leslie,  S.J.  Majumdar,  A.  Pang  and  F.  Lekien 
(2006)  Quantifying  Uncertainties  in  Ocean  Predictions.  Oceanography ,  Spe¬ 
cial  issue  on  “Advances  in  Computational  Oceanography” ,  T.  Paluszkiewicz 
and  S.  Harper,  Eds.,  19(1).  92-105. 

Lermusiaux,  P.F.J.  and  F.  Lekien  (2005)  Dynamics  and  Lagrangian  Co¬ 
herent  Structures  in  the  Ocean  and  their  Uncertainty.  Extended  Abstract 
in  Dynamical  System  Methods  in  Fluid  Dynamics ,  Oberwolfach  Workshop. 
J.E.  Marsden  and  J.  Scheurle  (Eds.),  Mathematisches  Forschungsinstitut 
Oberwolfach,  July  31st  -  August  6th,  2005,  Germany,  19-20. 

Lermusiaux  P.F.J. ,  A.R.  Robinson,  P.J.  Haley  and  W.G.  Leslie  (2002)  Ad¬ 
vanced  interdisciplinary  data  assimilation:  Filtering  and  smoothing  via  Er¬ 
ror  Subspace  Statistical  Estimation.  Proceedings  of  The  OCEANS  2002 
MTS/IEEE  conference,  Holland  Publications,  pp.  795-802. 

Li,  P.,  Y.  Chao,  Q.  Vu,  Z.  Li,  J.  Farrara,  H.  Zhang  and  X.  Wang  (2006) 
OurOcean  -  An  Integrated  Solution  to  Ocean  Monitoring  and  Forecasting. 
Proceedings  of  2006  Oceans  Conference ,  Boston,  September  2006. 

Liang,  X.S.  and  A.R.  Robinson  (2008)  Multiscale  Processes  and  Dynamics  for 
the  Monterey  Bay  Region  Circulation.  Journal  of  Physical  Oceanography. 
In  Press. 

Lozano,  C.J.,  P.J.  Haley,  H.G.  Arango,  Q.  Sloan  and  A.R.  Robinson  (1994) 
Harvard  coastal/deep  water  Primitive  Equation  model.  Harvard  Open 
Ocean  Model  Reports,  52,  Harvard  University,  Cambridge,  MA. 

Lozano,  C.J.,  A.R.  Robinson,  H.G.  Arango,  A.  Gangopadhyay,  N.Q.  Sloan, 
P.J.  Haley,  Jr.  and  W.G.  Leslie  (1996)  An  Interdisciplinary  Ocean 
Prediction  System:  Assimilation  Strategies  and  Structured  Data  Mod¬ 
els.  In  Modern  Approaches  to  Data  Assimilation  in  Ocean  Modelling  (P. 
Malanotte-Rizzoli,  Editors),  Elsevier  Oceanography  Series,  Elsevier  Science, 
The  Netherlands.  413-452. 

Lynn,  R.J.  and  J.J.  Simpson  (1987)  The  California  Current  system:  the  sea¬ 
sonal  variability  of  its  physical  characteristics.  J.  Geophys  Res.,  92.  12,947- 
12,966. 

Marchesicllo,  P.,  J.C.  McWilliams,  and  A.  Shchepetkin  (2003)  Equilib¬ 
rium  structure  and  dynamics  of  the  California  Current  System.  J.  Phys. 
Oceanogr.,  33.  753-783. 

MBARI  (2000)  Monterey  Bay  Multibeam  Survey;  Digital  Data  Series  No.  3. 
©2000  MBARI. 

Orlanski,  I.  (1976)  A  Simple  Boundary  Condition  for  Unbounded  Hyperbolic 
Flows.  J.  Comp.  Phys.,  21.  251-269. 

Pacanowski,  R.C.  and  G.H.  Philander  (1981)  Parameterization  of  Vertical 
Mixing  in  Numerical  Models  of  Tropical  Ocean.  J.  Phys.  Oceanogr .,  11, 
1443-1451. 

Pierce,  S.  D.,  R.L.  Smith,  P.M.  Kosro,  J.A.  Barth,  and  C.D.  Wilson  (2000) 
Continuity  of  the  poleward  undercurrent  along  the  eastern  boundary  of  the 
mid-latitude  Pacific.  Deep  Sea  Res.  II,  47(5-6).  811-829. 

Ramp,  S.R.,  R.E.  Davis,  N.E.  Leonard,  I.  Shulman,  Y.  Chao,  A.R.  Robinson, 


29 


J.  Marsden,  P.  Lermusiaux,  D.  Frantantoni,  J.D.  Paduan,  F.  Chavez,  X.S. 
Liang,  W.  Leslie  and  Z.  Li  (2008)  The  Autonomous  Ocean  Sampling 
Network  (AOSN)  Predictive  Skill  Experiment  in  the  Monterey  Bay.  Deep- 
Sea  Res.  II,  (submitted). 

Ramp,  S.R.,  P.F.  Jessen,  K.H.  Brink,  P.P.  Niiler,  F.L.  Daggett  and  J.S.  Best 
(1991)  The  physical  structure  of  cold  filaments  near  Point  Arena,  California, 
During  June  1987.  J.  Geophys.  Res.,  96.  14,859-14,883. 

Ramp,  S.R.,  J.D.  Paduan,  I.  Shulman,  J.  Kindle,  F.L.  Balir  and  F.  Chavez 
(2005)  Observations  of  Upwelling  and  Relaxation  Events  in  the  Northern 
Monterey  Bay  during  August  2000.  Journal  of  Geophysical  Research,  110, 
C07013. 

Robinson,  A.R.  (1999)  Forecasting  and  simulating  coastal  ocean  processes 
and  variabilities  with  the  Harvard  Ocean  Prediction  System.  In:  Coastal 
Ocean  Prediction,  (C.N.K.  Mooers,  ed.),  AGLT  Coastal  and  Estuarine  Stud¬ 
ies  Series.  American  Geophysical  Union.  77-100. 

Robinson,  A.R.,  H.G.  Arango,  A.  Warn-Varnas,  W.G.  Leslie,  A.J.  Miller,  P.J. 
Haley,  Jr.  and  C.J.  Lozano  (1996)  Real-Time  Regional  Forecasting.  Modern 
Approaches  to  Data  Assimilation  in  Ocean  Modelling  (P.  Malanotte-Rizzoli, 
Editors),  Elsevier  Oceanography  Series,  Elsevier  Science,  The  Netherlands. 
377-412. 

Robinson,  A.R.,  P.J.  Haley,  P.F.J.  Lermusiaux  and  W.G.  Leslie  (2002)  Pre¬ 
dictive  skill,  predictive  capability  and  predictability  in  ocean  forecasting.  In 
Proceedings  of  the  “The  OCEANS  2002  MTS/IEEE”  conference,  Holland 
Publications.  787-794. 

Robinson,  A.R.,  P.F.J.  Lermusiaux  and  N.Q.  Sloan  (1998)  Data  Assimilation. 
In  THE  SEA:  The  Global  Coastal  Ocean,  Volume  10:Processes  and  Methods 
(K.H.  Brink  and  A.R.  Robinson,  Editors),  John  Wiley  and  Sons,  NY,  541- 
594. 

Rosenfeld,  L.K.,  F.B.  Schwing,  N.  Garfield  and  D.E.  Tracy  (1994)  Bifur¬ 
cated  flow  from  an  upwelling  center:  a  cold  water  source  for  Monterey  Bay. 
Continental  Shelf  Res.,  14.  931-964. 

Rosenfeld,  L.,  I.  Shulman,  M.  Cook,  J.  Paduan,  and  L.  Shulman  (2008)  Eval¬ 
uation  of  a  regional  tidal  model  with  application  to  Central  California.  Deep 
Sea  Research  II,  (submitted). 

Rudnick,  D.L.,  R.E.  Davis,  C.C.  Eriksen,  D.M.  Fratantoni  and  M.J.  Perry 
(2004)  LInderwater  Gliders  for  Ocean  Research.  Journal  of  the  Marine 
Technology  Society.  38(2).  73-85. 

Scripps  Institution  of  Oceanography,  LIniversity  of  California  (1999)  Physical, 
chemical  and  biological  data,  CalCOFI  cruises  9802,  9803,  9804,  9805  and 
9806.  SIO  Ref.  99-9,  160pp. 

Shadden,  S.C.,  F.  Lekien  and  J.E.  Marsden  (2005)  Definition  and  properties 
of  Lagrangian  coherent  structures  from  finite-time  Lyapunov  exponents  in 
two-dimensional  aperiodic  flows.  Physica  D,  212(3-4).  271-304. 

Shadden,  S.C.,  F.  Lekien,  J.E.  Marsden  (2007)  Definition  and  Properties  of 
Lagrangian  Coherent  Structures:  Mixing  and  Transport  in  Two-Dimensional 


30 


Aperiodic  Flows.  Physica  D.  In  Press. 

Shadden,  S.C.,  F.  Lekien,  J.D.  Paduan,  F.P.  Chavez,  and  J.E.  Marsden  (2008) 
The  correlation  between  surface  drifters  and  coherent  structures  based  on 
high-frequency  radar  data  in  Monterey  Bay.  Deep  Sea  Research  II,  (submit¬ 
ted). 

Shapiro,  R.  (1970)  Smoothing,  Filtering  and  Boundary  Effects.  Reviews  of 
Geophysics  and  Space  Physics ,  8(2),  359-387. 

Sherman,  J.,  R.A.  Davis,  W.B.  Owens,  and  J.  Valdes  (2001)  The  autonomous 
underwater  glider  “Spray”,”  IEEE  Journal  of  Oceanic  Engineering ,  26(4). 
437-446,  Special  Issue  on  Autonomous  Ocean-Sampling  Networks. 

Shulman,  I.,  and  J.D.  Paduan  (2008)  Assimilation  of  HF  radar-derived  ra- 
dials  and  total  currents  in  the  Monterey  Bay  area.  Deep  Sea  Research  II, 
(submitted). 

Shulman,  I.,  C.  Rowley,  S.  Anderson,  S.  DeRada,  J.  Kindle,  P.  Martin,  J. 
Doyle,  J.  Cummings,  S.  Ramp,  F.  Chavez,  D.  Fratantoni,  and  R.  Davis 
(2008)  Impact  of  Glider  Data  Assimilation  on  Coastal  Model  Predictions, 
Deep  Sea  Research  II,  (submitted). 

Shulman,  I.,  C.R.  Wu,  J.K.  Lewis,  J.D.  Paduan,  L.K.  Rosenfeld,  J.C.  Kindle, 
S.R.  Ramp  and  C.A.  Collins  (2002)  High  resolution  Modeling  and  Data 
Assimilation  in  the  Monterey  Bay  Area,  Continental  Shelf  Research,  22, 
1129-1151. 

Sloan,  N.Q.,  III  (1996)  Dynamics  of  a  Shelf-Slope  Front:  Process  Studies 
and  Data-Driven  Simulations  in  the  Middle  Atlantic  Bight.  PhD  Thesis, 
Harvard  LIniversity,  Cambridge,  MA. 

Spall,  M.A.  and  W.R.  Holland  (1991)  A  Nested  Primitive  Equation  Model 
for  Oceanic  Applications.  J.  Phys.  Oceanog.,  21.  205-220. 

Spall,  M.A.  and  A.R.  Robinson  (1989)  A  New  Open  Ocean,  Hybrid  Coor¬ 
dinate  Primitive  Equation  Model.  Mathematics  and  Computers  in  Simula¬ 
tions,  31,  241-269. 

Steed,  C.A.,  J.E.  Braud,  and  K.A.  Koehler  (2002)  VGRID:  a  generic,  dynamic 
HDF5  storage  model  for  georeferenced  grid  data.  Oceans  ’02  MTS/IEEE, 
2.  900-907. 

Strub,  P.T.  and  C.  James  (2000)  Altimeter-derived  variability  of  surface 
velocities  in  the  California  Current  System:  2.  Seasonal  circulation  and  eddy 
statistics.  Deep-Sea  Res.  II,  45.  831-870. 

Strub,  P.T.,  P.M.  Kosro  and  A.  Huyer  (1991)  The  nature  of  the  cold  filaments 
in  the  California  Current  system.  J.  Geophys.  Res.,  96.  14,743-14,768. 

Traganza,  E.D.,  J.C.  Conrad,  L.C.  Breaker,  and  F.A.  Richards  (1981)  Satel¬ 
lite  observations  of  a  cyclonic  upwelling  system  and  giant  plume  in  the 
California  Current.  In  Coastal  Upwelling,  DOE  Symposium  on  Coastal  Up- 
welling,  AGU.  228-241 

Tukey,  J.W.  (1977)  Exploratory  Data  Analysis.  Addison- Wesley,  Reading 
MA. 

Wang,  X.,  Y.  Chao,  C.  Dong,  J.  Farrara,  Z.  Li,  J.C.  McWilliams,  J.D.  Paduan, 
and  L.K.  Rosenfeld  (2008)  Modeling  Tides  in  Monterey  Bay,  California. 


31 


Deep  Sea  Research  II,  (submitted). 

Wang,  P.,  Y.T.  Song,  Y.  Chao,  H.  Zhang  (2005)  Parallel  Computation  of 
the  Regional  Ocean  Modeling  System,  The  International  Journal  of  High 
Performance  Computing  Applications,  19(4).  375-385. 

Wickham,  J.B.  (1975)  Observations  of  the  California  countercurrent.  J.  Mar. 
Res.,  33.  325-340. 


32 


List  of  Figures 


1  Modeling  domains  for  AOSN-II  with  bathymetry.  Top  panel 
shows  the  “Data  Domain”  inside  of  the  larger  “Offshore 
Domain” .  Bottom  panel  shows  a  zoom  of  the  “Data  Domain” . 

Also  shown  are  the  main  dynamical  features:  upwelling  centers 
at  Point  Ano  Nuevo  (AN)  and  Point  Sur  (PS)  (blue);  coastal 
current,  eddies,  filaments,  etc.  (black);  California  Undercurrent 
(CUC)  (green);  California  Current  (CC)  (magenta).  36 

2  Wind  vectors  (knots)  at  the  Ml  and  M2  moorings.  37 

3  Positions  of  primary  sources  of  in  situ  data,  (a-c):  locations  of 

the  R/V  Pt.  Sur  (repeat)  CTD  surveys,  (d):  locations  of  the 
pseudo-profiles  from  the  Woods  Hole  and  Scripps  gliders.  38 

4  Flowchart  of  the  canonical  day’s  operations  for  AOSN-II.  The 
tuning  operations  (shown  in  dashed  lined)  were  not  performed 
in  AOSN-II  but  were  suggested  by  the  reanalysis.  These 
tunings  are  a  form  of  adaptive  modeling  (Lermusiaux,  2007) 
in  which  a  series  of  sensitivity  studies  are  used  to  improve 
the  parameterizations  during  the  operation.  This  notion  was 
successfully  employed  during  the  subsequent  Monterey  Bay 

2006  experiment.  39 

5  Temperature  (°C)  and  velocity  (cm/s)  fields  for  Aug  8,  2003 

issued  during  the  real-time  exercise.  These  fields  show  the 
general  structures  that  were  observed  in  real-time.  40 

6  Comparison  of  temperature  (°C)  and  velocity  (cm/s)  for 

stand-alone,  nested  and  ESSE  forecasts  for  13  Aug  2003.  41 

7  Comparisons  of  (top  row)  available  real-time  SST  (°C) 

images  from  satellite  (Aug  12  and  21)  and  aircraft  (Aug  13) 
with  (bottom  row)  corresponding  real-time  forecast  surface 
temperature  (°C)  fields.  Note  timing  conventions  (e.g.  0000Z 
13  Aug  forecast  corresponds  to  1600  local  12  Aug  and  is 
compared  to  12  Aug  image).  42 

8  Comparisons  of  (top  row)  available  real-time  CODAR  (cm/s) 

images  with  (bottom  row)  corresponding  real-time  forecast 
surface  velocity  (cm/s)  fields.  43 


33 


9  Real-time  volume-averaged  RMS  error  differences.  Persistence 
RMS  -  Forecast  RMS,  a  positive  result  indicates  skill.  The 
left  column  shows  results  for  temperature  (°C),  the  right  for 
salinity  (PSU).  The  top  row  shows  results  for  1  day  forecasts, 

the  bottom  row  for  2  day  forecasts.  44 

10  Real-time  Temperature  RMS  differences  (°C)  for  1  day 

forecasts.  Persistence  RMS  -  Forecast  RMS,  a  positive  result 
indicates  skill.  45 

11  Real-time  Temperature  RMS  differences  (°C)  for  2  day 

forecasts.  Persistence  RMS  -  Forecast  RMS,  a  positive  result 
indicates  skill.  46 

12  Real-time  volume-averaged  PCC  differences.  Forecast  PCC 
-  Persistence  PCC,  a  positive  result  indicates  skill.  The  left 
column  shows  results  for  temperature,  the  right  for  salinity. 

The  top  row  shows  results  for  1  day  forecasts,  the  bottom  row 

for  2  day  forecasts.  47 

13  Reanalysis  surface  temperature  (°C)  and  velocity  (cm/s)  every 

2  days  for  the  period  10  Aug  -  1  Sep,  2003.  48 

14  Reanalysis  30m  temperature  (°C)  and  velocity  (cm/s)  every  2 

days  for  the  period  10  Aug  -  1  Sep,  2003.  49 

15  Reanalysis  Temperature  RMS  differences  (°C).  Persistence 
RMS  -  Reanalysis  RMS,  a  positive  result  indicates  skill.  In  this 
reanalysis  simulation,  only  data  up  to  Aug  20  was  assimilated. 
Persistence  is  defined  as  the  Aug  20  objective  analysis  of  all 

data  up  to  Aug  20.  50 

16  Reanalysis  Temperature  Bias  (°C)  and  RMS  (°C)  comparisons, 
averaged  over  Aug  21-25.  In  this  reanalysis  simulation,  only 
data  up  to  Aug  20  was  assimilated.  Persistence  is  defined  as 

the  Aug  20  objective  analysis  of  all  data  up  to  Aug  20.  51 
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Initialization 

Assimilation 

Synoptic 

Mean 

Synoptic 

Mean 

Decay  (km) 

15 

45 

5 

25 

Zero-Cross  (km) 

37.5 

112.5 

12.5 

50 

Time-Decay  (day) 

10 

1000 

1 

80 

Table  1 

Objective  Analysis  correlation  parameters. 


Time 

Weights  (T,S) 

Weights  (Uinternal) 

tassim  ~  0.25 

0.333 

0.167 

tassim  ~  0.125 

0.666 

0.333 

tassim 

0.999 

0.5 

tassim  0.25 

0.333 

0.167 

Table  2 

Optimal  Interpolation  assimilation  ramping  parameters. 


Boundary  k  fix 

T\ 

(days) 

T2 

(days) 

7 

Real-Time  Values 

w 

5 

1 

4 

5 

s 

5 

0.1 

1.5 

5 

n 

5 

0.1 

3.9 

6 

Reanalysis  Values 

w 

5 

1 

4 

5 

s 

5 

1 

3.9 

5 

n 

5 

1 

3.9 

6 

Table  3 

AOSN-BC  parameters  used  in  real-time  forecasts  and  reanalysis. 
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Fig.  1.  Modeling  domains  for  AOSN-II  with  bathymetry.  Top  panel  shows  the  “Data 
Domain”  inside  of  the  larger  “Offshore  Domain”.  Bottom  panel  shows  a  zoom  of 
the  “Data  Domain”.  Also  shown  are  the  main  dynamical  features:  upwelling  cen¬ 
ters  at  Point  Aho  Nuevo  (AN)  and  Point  Sur  (PS)  (blue);  coastal  current,  eddies, 
filaments,  etc.  (black);  California  Undercurrent  (CUC)  (green);  California  Current 
(CC)  (magenta). 
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08/03 


Wind  Vectors  at  Buoy  Ml 
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Fig.  2.  Wind  vectors  (knots)  at  the  Ml  and  M2  moorings. 
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(a) 


Point  Sur  Data:  2-6  August  2003 


(c) 

Point  Sur  Data:  3-6  September  2003 


(b) 


Point  Sur  Data:  21-25  August  2003 


(d) 

AOSN-II  Glider  Positions 


-123  -122.5  -122  -121.5 


Longitude 


Fig.  3.  Positions  of  primary  sources  of  in  situ  data,  (a-c):  locations  of  the  R/V  Pt. 
Sur  (repeat)  CTD  surveys,  (d):  locations  of  the  pseudo-profiles  from  the  Woods 
Hole  and  Scripps  gliders. 
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Fig.  4.  Flowchart  of  the  canonical  day’s  operations  for  AOSN-II.  The  tuning  opera¬ 
tions  (shown  in  dashed  lined)  were  not  performed  in  AOSN-II  but  were  suggested  by 
the  reanalysis.  These  tunings  are  a  form  of  adaptive  modeling  (Lermusiaux,  2007) 
in  which  a  series  of  sensitivity  studies  are  used  to  improve  the  parameterizations 
during  the  operation.  This  notion  was  successfully  employed  during  the  subsequent 
Monterey  Bay  2006  experiment. 
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Fig.  5.  Temperature  (°C)  and  velocity  (cm/s)  fields  for  Aug  8,  2003  issued  during 
the  real-time  exercise.  These  fields  show  the  general  structures  that  were  observed 
in  real-time. 


40 


Stand  Alone  OI  Nested  OI  Stand  Alone  ESSE 


0.00.00 Z  '3  Aug  2003 
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Fig.  6.  Comparison  of  temperature  (°C)  and  velocity  (cnr/s)  for  stand-alone,  nested 
and  ESSE  forecasts  for  13  Aug  2003. 
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Fig.  7.  Comparisons  of  (top  row)  available  real-time  SST  (°C)  images  from  satellite  (Aug  12  and  21)  and  aircraft  (Aug  13)  with  (bottom 
row)  corresponding  real-time  forecast  surface  temperature  (°C)  fields.  Note  timing  conventions  (e.g.  0000Z  13  Aug  forecast  corresponds 
to  1600  local  12  Aug  and  is  compared  to  12  Aug  image). 
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Fig.  8.  Comparisons  of  (top  row)  available  real-time  CODAR  (cm/s)  images  with  (bottom  row)  corresponding  real-time  forecast  surface 
velocity  (cm/s)  fields. 
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Fig.  9.  Real-time  volume-averaged  RMS  error  differences.  Persistence  RMS  -  Fore¬ 
cast  RMS,  a  positive  result  indicates  skill.  The  left  column  shows  results  for  tem¬ 
perature  (°C),  the  right  for  salinity  (PSU).  The  top  row  shows  results  for  1  day 
forecasts,  the  bottom  row  for  2  day  forecasts. 
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Fig.  10.  Real-time  Temperature  RMS  differences  (°C)  for  1  day  forecasts.  Persistence 
RMS  -  Forecast  RMS,  a  positive  result  indicates  skill. 
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Fig.  11.  Real-time  Temperature  RMS  differences  (°C)  for  2  day  forecasts.  Persistence 
RMS  -  Forecast  RMS,  a  positive  result  indicates  skill. 
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Fig.  12.  Real-time  volume-averaged  PCC  differences.  Forecast  PCC  -  Persistence 
PCC,  a  positive  result  indicates  skill.  The  left  column  shows  results  for  temperature, 
the  right  for  salinity.  The  top  row  shows  results  for  1  day  forecasts,  the  bottom  row 
for  2  day  forecasts. 
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Fig.  13.  Reanalysis  surface  temperature  (°C)  and  velocity  (cm/s)  every  2  days  for 
the  period  10  Aug  -  1  Sep,  2003. 
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Fig.  14.  Reanalysis  30m  temperature  (°C)  and  velocity  (cm/s)  every  2  days  for  the 
period  10  Aug  -  1  Sep,  2003. 


49 


Trms 

Persistence 


-  Trms„ 


Reanalysis 


08/21  08/22  08/23  08/24  08/25 


0 


08/21  08/22  08/23  08/24  08/25 


1.5 


08/21 

08/22 

08/23 

08/24 

08/25 

1.5 

- 

1 

- 

0  5 

_ 

- © - — 

0 

08/21 

08/22 

08/23 

08/24 

08/25 

1.5 


08/21  08/22  08/23  08/24  08/25 


Fig.  15.  Reanalysis  Temperature  RMS  differences  (°C).  Persistence  RMS  -  Reanaly¬ 
sis  RMS,  a  positive  result  indicates  skill.  In  this  reanalysis  simulation,  only  data  up 
to  Aug  20  was  assimilated.  Persistence  is  defined  as  the  Aug  20  objective  analysis 
of  all  data  up  to  Aug  20. 
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Fig.  16.  Reanalysis  Temperature  Bias  (°C)  and  RMS  (°C)  comparisons,  averaged 
over  Aug  21-25.  In  this  reanalysis  simulation,  only  data  up  to  Aug  20  was  assimi¬ 
lated.  Persistence  is  defined  as  the  Aug  20  objective  analysis  of  all  data  up  to  Aug 
20. 
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